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While, to ensure successful cosmology, dark matter (DM) must kinematically decouple from the 
standard model plasma very early in the history of the Universe, it can remain coupled to a bath 
of "dark radiation" until a relatively late epoch. One minimal theory that realizes such a scenario 
is the Atomic Dark Matter model, in which two fermions oppositely charged under a new U(l) 
dark force are initially coupled to a thermal bath of "dark photons" but eventually recombine 
into neutral atom-like bound states and begin forming gravitationally-bound structures. As dark 
atoms have (dark) atom-sized geometric cross sections, this model also provides an example of 
self-interacting DM with a velocity-dependent cross section. Delayed kinetic decoupling in this 
scenario predicts novel DM properties on small scales but retains the success of cold DM on larger 
scales. We calculate the atomic physics necessary to capture the thermal history of this dark sector 
and show significant improvements over the standard atomic hydrogen calculation are needed. We 
solve the Boltzmann equations that govern the evolution of cosmological fluctuations in this model 
and find in detail the impact of the atomic DM scenario on the matter power spectrum and the 
cosmic microwave background (CMB). This scenario imprints a new length scale, the Dark- Acoustic- 
Oscillation (DAO) scale, on the matter density field. This DAO scale shapes the small-scale matter 
power spectrum and determines the minimal DM halo mass at late times which may be many orders 
of magnitude larger than in a typical WIMP scenario. This model necessarily includes an extra dark 
radiation component, which may be favoured by current CMB experiments, and we quantify CMB 
signatures that distinguish an atomic DM scenario from a standard ACDM model containing extra 
free-streaming particles. We finally discuss the impacts of atomic DM on galactic dynamics and 
show that these provide the strongest constraints on the model. 

PACS numbers: 98.80.-k,98.80. Jk 



I. INTRODUCTION 

Although dark matter (DM) has been known to exist 
for several decades [1 - .-3 J , its physical nature remains one 
of the deepest mysteries of modern science. Observations 
show that DM is mostly cold, collisionless and interacts 
very weakly (if at all) with standard-model (SM) parti- 
cles. Many models that fit this picture very well have 
been proposed through the years, including the promi- 
nent weakly-interacting-massive-particle (WIMP) mod- 
els. Distinguishing between these different scenarios is 
crucial if we have any hope to pinpoint the nature of 
DM. In this respect, models incorporating new physics 
that predict novel observational signatures have a clear 
advantage in disentangling the DM puzzle. 

While the cold dark-matter (CDM) paradigm [U [5] 
has been extremely successful at describing observations 
from cosmological scales to galactic scales, recent ob- 
servations of small nearby galaxies seem to be in ten- 
sion with the CDM scenario. In addition to the so- 
called "missing-satellite problem" [3H5] which refers to 
the apparent under-abundance of light Milky- Way satel- 
lites, observations of the inner mass profiles of Dwarf 
Spheroidal (dSph) galaxies indicate that they are consis- 
tent with a core while CDM simulations favor a cuspier 
profile [9lU3|. Further, it has been pointed out recently 
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that the most massive galactic subhalos from CDM sim- 
ulations are too dense to host the brightest Milky Way 
satellites [TU [15] . While it is plausible that these prob- 
lems could be alleviated by including the appropriate 
baryonic physics in CDM simulations [TBT - 12?] . these ob- 
servations may be pointing toward physics beyond the 
vanilla CDM paradigm. 

We adopt here the point of view that the above ten- 
sions between the ACDM paradigm and the astrophysical 
observations may be resolved by modifying the micro- 
physics of the DM sector. Various scenarios have been 
proposed along those lines in the literature, most of which 
could either be classified as interacting DM models [28] . 
warm DM scenarios [2"M3"!2] , or even hot DM models (see 
e.g. [33]). The former is tightly constrained by the ob- 
served ellipticity of DM halos [SHE!] and by the apparent 
survivability of DM halos to evaporation in clusters [33] , 
while the latter is unlikely to be able to address all of the 
CDM issues discussed above [37] . It was realized recently 
[3^H1U] that interacting DM with a velocity-dependent 
cross section could avoid the ellipticity and evaporation 
constraints while alleviating the tension between the sim- 
ulations and the Milky- Way satellites. 

In this paper, we investigate in detail the cosmology of 
a dark sector made of hidden hydrogen-like bound states 
[3TH33] . This so-called "Atomic Dark Matter" model re- 
tains the success of CDM on large cosmological scales but 
modifies the DM dynamics on sub-galactic scales. Sim- 
ilar to conventional atoms, atomic DM is kinematically 
coupled to a thermal bath of "dark" radiation (DR) until 
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late times, which significantly delays the growth of mat- 
ter perturbations on small scales. During the decoupling 
epoch, diffusion and acoustic damping substantially re- 
duce the amplitude of sub-horizon perturbations, effec- 
tively wiping out structures on scales below this damp- 
ing horizon, hence mimicking the effect of free-streaming. 
This model thus naturally provides a way to suppress the 
faint end of the galaxy luminosity function. After kine- 
matic decoupling, the acoustic oscillations of the dark 
plasma remain imprinted on the small-scale matter power 
spectrum. Note that these physical processes also take 
place in a canonical WIMP scenario [45 47J . The crucial 
difference here is that the kinematic decoupling of atomic 
DM typically happens much later than that of a WIMP 
and therefore can have an impact on much larger, if not 
observable, scales. 

Beyond its effect on the matter power spectrum, the 
atomic DM scenario also impacts the cosmic microwave 
background (CMB) through the effects of the DR on 
the amplitude and phases of cosmological perturbations. 
We identify a key CMB signature that distinguishes the 
atomic DM scenario from a simple ACDM model incorpo- 
rating extra relativistic degrees of freedom. Indeed, while 
models incorporating extra radiation always assume it 
to free-stream like neutrinos, the DR in the atomic DM 
model can only start free-streaming after it decouples 
from the DM. Cosmological modes entering the Hubble 
horizon before the onset of the free-streaming display a 
slightly different behavior than Fourier modes crossing 
the causal horizon after the decoupling of the DR. 

Since dark atoms have a much larger geometric cross 
section than point particles, this scenario falls into 
the category of interacting DM models. Interestingly, 
the atomic physics naturally gives rise to a velocity- 
dependent interaction cross section. It is therefore pos- 
sible that atomic DM could address some of the dSph- 
galaxy problems while evading the ellipticity and evap- 
oration constraints. In the following, we focus on the 
ellipticity constraint and show that while it strongly con- 
strains the parameter space of the model, there are plenty 
of parameter values for which the galactic dynamics is 
unaffected. 

For clarity and completeness, it is important to men- 
tion that the term "dark atoms" has been used in various 
contexts in the literature. For example, dark atoms nat- 
urally appears in Mirror DM models [48-452"] . In [53] |53] , 
the term "dark atom" refers to a bound state between 
a new stable particle charged under ordinary electro- 
magnetism and a helium nucleus. The authors of [55] 
explored a model in which dark matter is made of su- 
persymmetric dark atoms. In the present work, the 
phrase "dark atom" is exclusively used to describe the 
bound state of two dark fermions (i.e. neutral under the 
SM gauge group) oppositely charged under a new gauge 
U{\)d interaction. 

This paper is organized as follow. In section [TTJ we 
introduce the atomic DM model and discuss the param- 
eters necessary to describe the theory. In section |III| we 



analyze the thermal history of atomic DM from its ini- 
tial hot plasma state to its cold and collisionless stage 
at late times. We include a thorough discussion of the 
thermal and kinetic decoupling epoch and present an in- 
depth analysis of the dark-recombination process. In sec- 
tion |IV] we solve the modified cosmological perturbation 
equations and discuss the different regimes that a DM 
fluctuation encounters. In section [V] we present an anal- 
ysis of the new features in the matter power spectrum 
and in the CMB due to the dark atomic physics. We also 
discuss the Lyman-a (Ly-a) forest bounds on the pa- 
rameter space of atomic DM. In section [yl] we consider 
the stringent astrophysical constraints on the model that 
are obtained by imposing that DM is effectively collision- 
less inside galactic halos. We revisit in section |VII| the 
direct-detection signatures of atomic DM proposed in the 
literature in light of our new analysis. We finally discuss 
our results in section |VIII| and point out possible new 
avenues of research. 



II. THE MODEL 

The atomic DM model [HTH5] is composed of two 
oppositely-charged massive fermions interacting through 
a dark massless (or nearly massless) U{X)d gauge bo- 
sorfj] In analogy with the regular baryonic sector, we call 
the lighter fermion "dark electron" , the heavier fermion 
"dark proton" and the massless gauge boson "dark pho- 
ton". In this work, we assume that the DM relic abun- 
dance is set by some UV physics in the early Universe. 
For example, the dark sector could contain an asymmetry 
that fixes the DM relic density (see e.g. [55]). Since we 
are mostly interested in the late-time cosmological and 
astrophysical impacts of atomic DM, we do not expect 
the details of this high-energy completion to play any 
role in our results. See [H] 1571450] for examples of DM 
production mechanism that can lead to atom-like DM 
particles. 

Four parameters are necessary to fully characterize the 
physics of the dark sector in this model. These are the 
DM mass (i.e. the mass of the bound state) mo, the dark 
fine-structure constant ao, the binding energy of the 
bound state Bjj and the present-day ratio of the DR tem- 
perature to the CMB temperature £ = (Td/Tcmb)U=o- 
Other combinations of parameters are possible but this 
particular set is physically transparent since Bd/£, fixes 
the redshift of dark recombination, mjy fixes the number- 
density of DM particles and a d governs the microscopic 
interactions between the dark sector constituents. These 



We note that dark atoms do not have to be made of two spin- 1/2 
fermions. For instance, they could arise from a bound state of 
two scalar particles or from a bound state of a spin-1/2 fermion 
with a scalar particle. 



3 



parameters are subject to the consistency constraint 



which ensures that m e p < (mrj + -Bd)/2. This bound 
is saturated when the two fermions have equal masses. 
Here, m e and m p stand for the dark-electron mass and 
the dark-proton mass, respectively. We give their values 
in terms of mo, Bo, and an in appendix |A"} 

On a more general level, the atomic DM scenario can 
be considered as a toy model of a more complete the- 
ory involving a hidden dark plasma in the early Uni- 
verse. Indeed, atomic DM contains the all key ingredients 
of a dark-plasma theory (dark radiation, multiple par- 
ticles, kinetic and thermal decoupling, modified growth 
of DM fluctuations, long-range and short-range interac- 
tions, etc) with only minimal physical inputs. As such, 
the results presented in this paper should be understood 
in the broader context of a generalized dark-plasma the- 
ory ( see e.g. [61] for general cosmological constraints on 
this type of model) . 

Interestingly, the atomic DM scenario naturally en- 
globes the hidden charged DM models discussed in 
[6"2l 16"3"] as special cases. Moreover, in the limit of very 
large atomic binding energy and large dark fine-structure 
constant, dark atoms become basically undistinguishable 
from standard CDM particles. Therefore, the atomic DM 
scenario is a rather general testbed for physics beyond the 
vanilla WIMP/CDM paradigm. 



III. THERMAL HISTORY 

In the early Universe, the dark sector forms a tightly- 
coupled plasma much like the standard baryon-photon 
plasma. As the Universe cools down, three important 
transitions need to carefully be taken into account. First, 
once the dark sector temperature falls below Bp, it be- 
comes energetically favorable for the dark fermions to 
recombine into neutral dark atoms. Second, once the 
momentum transfer rate between the DR and the dark 
fermions falls below the Hubble expansion rate (kinetic 
decoupling), the DM effectively ceases to be dragged 
along by the radiation and can start to clump and form 
structures. Finally, once the energy transfer rate between 
the DM and the DR falls below the Hubble rate (thermal 
decoupling), the DM temperature Tjjm ceases to track 
that of DR and start cooling adiabatically. Accurately 
capturing these transitions and computing their impact 
on cosmological observables is a major goal of this paper. 

We begin this section by determining the big-bang nu- 
cleosynthesis (BBN) bound on the dark-photon temper- 
ature. We then discuss the recombination of dark atoms 
and their thermal coupling to the DR, emphasizing the 
differences between dark atoms and regular atomic hy- 
drogen. We finally present the solutions to the joint evo- 
lution of the dark-atom ionized fraction and temperature. 




FIG. 1. Effective number of dark sector relativistic degrees of 
freedom at the time of nucleosynthesis as a function of ob and 
Bd for dark atoms with mass mo = 1 GeV. Here, we have 
fixed £bbn = 0.5. We also display the consistency constraint 
given by Eq.[l] above which dark atoms do not exist. 

A. BBN Limit on Dark-Sector Temperature 

Observations of the relative abundance of light ele- 
ments put a bound on the possible number of relativistic 
degrees of freedom at the time of BBN. This limit is 
usually quoted in terms of the effective number of light 
neutrino species in thermal equilibrium at BBN; here we 
shall use the conservative estimate N u = 3.24 ± 1.2 (95% 
confidence) derived in Ref. [M]. More recent estimates 
[6"S"H6"9"] of N u are statistically consistent with this value. 
Assuming that the dark sector contributes gf D N relativis- 
tic degrees of freedom during BBN and further assuming 
three species of SM neutrinos, we obtain the bound 

9™d&bs < 2-52, (2) 

where £bbn is the ratio of the dark sector and visible 
temperatures at the time of nucleosynthesis. In the min- 
imal atomic DM scenario considered in this work, DM 
is totally decoupled from SM particles and therefore the 
entropy of the dark sector and the visible sector are sep- 
arately conserved 

BBN f 3 „BBN 
y*S,USBBN _ y*S,vis /„\ 
today £3 today ' s ' 

9*S,D S 5*S,vis 

where g*s.D is the present-day effective number of de- 
grees of freedom contributing to the entropy of the dark 
sector, sd oc g^sM^D- For the simplest model of dark 
atoms considered here, we have g^sD = 2 (i.e. only 
dark photons contribute). Similarly, g*s,vis is the effec- 
tive number of degrees of freedom contributing to the 
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entropy of the visible sector. For the particle content 



of the SM, we have g 



BBN 
*S,vis 



10.75 and g 



today 
*S,vis 



3.94. 



During BBN, both dark photons and dark electrons (to- 
gether with their antiparticles) can contribute to .gfj 3 ^ 
(we assume that the dark proton is massive enough to be 
non-relativistic at the time of BBN) . These dark compo- 
nents are kept in thermal equilibrium through Compton 
scattering and we therefore always have gfg^ = gf f) N . 
Figure nl shows the dependence of gfgp on au and 
Bjj for aark atoms with mass = 1 GeV and for 
Cbbn = 0.5. We see that there is a large parameter space 
for which dark electrons are relativistic at BBN, leading 
to gfg^, — 11/2 for these models. At late times, the ra- 
tio of the dark sector and visible-sector temperatures is 
given by 



' today BBN \ ^'^ 
9*S,vis9*S,D \ 

BBN today I ' 



„BBN ■ 



(4) 



Note that if gfg^, = 11/2, that is, if the dark electrons 
and anti-electrons annihilate after BBN, then £ s» £bbn 
since both the visible and dark sector are reheated by the 
same amount in this case (assuming the dark electron is a 
Dirac fermion) . Substituting Eq. [4] into Eq. [2] we obtain 
the constraint 



(„BBNU/12 

\y*,D > 



<0.71 (95% confidence). 



(5) 



We display this constraint in Fig. [2] where we observe 
that £ > 1 is disfavored by at least 4 standard devia- 
tions if we consider BBN alone. Note that atomic DM 
models generally predict a different number of effective 
relativistic degrees of freedom at BBN than at the time of 
hydrogen recombination. Given a choice of dark param- 
eters (an, Bd, friD, £), we can compute 5*,d(Td) using 
Eq. Bl The evolution of the dark-radiation tempera- 



ture is then given by the implicit equation 



T D (z)=T t D od ^(l + z) 



' today 
9*,D 

g*,D(TD 



1/3 



(6) 



This equation can easily be solved iteratively by substi- 
tuting the zeroth order solution Td(z) = x^ day (l + z) 
into gsr^DiTjj). In practice however, the annihilation of 
dark electron and dark positron has very little effect on 
the cosmological evolution and to a very good approxi- 
mation, we can take £ to be constant. 



B. Dark Recombination 

Once T D -C B D , it becomes energetically favorable for 
the dark plasma to recombine into neutral dark atoms. 
Letting n e and rip denote the number densities of dark 
electrons and of dark protons (both free and those bound 
in dark atoms) respectively, we can define the ionization 
fraction of the dark plasma as = n e /no. 
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FIG. 2. Joint BBN constraints on the present-day dark sec- 
tor temperature and on the effective number of dark sector 
relativistic degrees of freedom at the time of nucleosynthesis. 
As indicated, we display contours corresponding to 1-, 2-, and 
3-cr constraints. 



The phenomenology of dark-atom recombination can 
significantly differ from that of standard atomic hydro- 
gen. In the latter case, atoms cannot directly recombine 
to their ground state (Case A) as this results in the emis- 
sion of a Lyman-continuum photon that can immediately 
ionize a neutral atom [70U71j . In the case of dark atoms, 
it can be shown that for all but the most weakly-coupled 
models, the dark sector is always optically thick to dark 
Lyman-continuum photons and dark atoms cannot effec- 
tively recombine directly to their ground state. Dark re- 
combination thus needs to proceed through excited states 
(Case B) for most models. For the very weakly-coupled 
models, the dark sector can become optically thin to the 
Lyman-continuum dark photons, hence allowing dark re- 
combination directly to the ground state. However, the 
very small coupling between the different dark sector con- 
stituents implies that very little recombination actually 
takes place for these models. We therefore conclude that 
whenever dark recombination is at all appreciable, it pro- 
ceeds through case B. Hence, we neglect ground-state re- 
combination in the following. 

We can classify dark-atom models in three broad cate- 
gories according to their recombination phenomenology. 
The first category encompasses models for which the ma- 
jority of the recombination process takes place in ther- 
modynamic equilibrium. This category includes models 
having relatively large values of the dark fine-structure 
constant (od > 0.1) and low-mass models for which the 
number density is very large during dark recombination. 
For these scenarios, the recombination history is very well 
captured by the Saha ionization equation until the re- 
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combination timescale becomes comparable to the Hub- 
ble expansion rate and the ionized fraction freezes out. 
In this case, the recombination process is insensitive to 
the exact population of the excited states and to the de- 
tails of the radiative transfer between the various atomic 
lines, making the determination of the ionization history 
quite simple. 

The second category encompasses models for the re- 
combination timescale is slightly shorter or comparable 
to the Hubble expansion rate. Standard atomic hydro- 
gen falls into this category. For these models, the details 
of the transitions between the excited states and ground 
state are relevant. The basic picture is as follow. Once 
a dark electron has reached an excited state, it rapidly 
cascades down to the n = 2 state. Dark electrons can 
then reach the ground state via the 2s-ls two-photon 
transition or by emitting a Lyman-a dark photon which 
redshifts out of the line wing. The recombination pro- 
cess is thus governed by the rate of escape of the dark 
Lyman-a photons out of their resonance line and by the 
rate of two-photon decay [70- 72 . While these two rates 
are comparable for regular atomic hydrogen, they may 
differ significantly in the case of dark atoms. In addi- 
tion, it has been shown [73HH3] that the details of the 
radiative transfer process can affect recombination at the 
few percent level for this type of models. While we do 
not take these corrections into account in this work, we 
expect them to have a negligible impact on our results. 

The last category englobes models for which the re- 
combination timescale is longer than the Hubble rate. 
Atomic DM scenarios having very small values of the 
dark fine-structure constant or very large masses fall into 
this category. Not surprisingly, the dark sector remains 
mostly ionized for these models. The small amount of 
neutral dark atoms that do form in this case is mostly 
controlled by the small value of the effective recombina- 
tion rate and is, to a good approximation, insensitive to 
the radiative transfer between the various atomic lines. 

For all models, we solve for the ionization history of 
the dark sector using the Effective Multi-Level Atom 
(EMLA) method presented in Refs. [55J US]- For our 
purpose, we consider a four-level effective dark atom con- 
sisting of the two interface states 2s and 2p, the ground 
state (Is), and a continuum. The recombination process 
is then governed by a set of rate equations of the form 



dx D _ x 2 D n D (A 2 £ 
dz 



At) 



B 2 D s x 2s 



B 2 £x 2p 



H{z){l + z) 



(7) 



dxo 



B 2 D s x 2s 



c 2 D n D A 



D 



n 



2p,2s 



(3x 2s - x 2p ) 



A 2 



,{x 2s 



H(z)(l + z) 
- {l-x D )e^ 



H(z)(l 



z 



(8) 



dx 2p _ B 2 ^x 2p - x 2 D n D A 2 £ + H 2p . 2s (x 2p ~ 3x 2s ) 
dz ~~ H(z)(l + z) 

Rhya(x 2p - 3(1 - x D )e it d ) 



H(z)(l + z) 



(9) 



where A^ = A^(Tdm,Td) is the effective recombina- 
tion coefficient to the interface state nl, Bp — B^(Td) 
is the effective photoionization rate of the interface state 
nl, H(z) is the Hubble expansion rate and z is the red- 
shift which we use here as a time variable. Ri_, ya is the 
rate at which dark Lyman-a photons escape the Lya res- 
onance due to the expansion of the Universe [S3] 



R 



Lya 



2 9 a 3 D B D (l-e-^) 



3 8 



TS 



(10) 



where rs is the Sobolev optical depth for the Lyman-a 
transition. Here and in the remainder of this paper, we 
set h — c — fee = 1 ■ The second factor in Eq. [10] is the 
probability that a Lyman-a dark photon redshifts out 
of the line wing before being reabsorbed. The Lyman-a 
Sobolev optical depth is given by [84] 



ts 



2 9 a 3 D B D n D (l - x D ) 



x 2p 



3(i-a;i>; 



(ii) 



where Vhya is the frequency of the Lyman-a dark pho- 
tons. A2 Sj is denotes the rate of the forbidden two-photon 
transition between the 2s state and the ground state. It 
is given by 



A 



2s- Is 



Bd 



3.22458 s" 



(12) 



Here a em is the SM fine-structure constant and -Bh — 
13.6 eV is the binding energy of regular hydrogen. lZ 2p . 2s 
is the effective transition rate between the interface states 
2 s and 2p and is given by [55] 



2p.2s 



2/>- 



l.Pn. 



(13) 



where R 2p -+ n i is the bound-bound transition rate be- 
tween the 2p state and the "interior" state nl (n > 2) 
and P n i^, 2s stands for the probability that a dark atom 
in a state nl will decay to the 2s state. Details of the 
computation of R 2p ^ n i and P„i^2s can be found in [85] . 

Previous studies of atomic DM have used an approxi- 
mate recombination coefficient given by [SOI HI] 



A D {T DM ) = 0.448 



64tt 



- 

/27^n 2 D \TdmJ 



1/2 



\Tdm) 

(14) 

This rate is problematic for three reasons. First, since 
it is only a function of the DM temperature, it fails to 
take into account the contribution from stimulated re- 
combination which explicitly depends on Tjj. Second, it 
does not take into account the effect of transitions among 
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the high-n atomic states and assume equilibrium among 
the angular momentum substates (an approximately 14% 
correction). Finally, it is inaccurate for Tdm > Bp and 
Tdm Bjj. To improve this situation and obtain an 
accurate picture of dark recombination, we compute a 
new recombination coefficient using the method outlined 
in |85j . Our recombination coefficient is given by 



DM, Td 



Oinl 



^ On'l'Pn 



(15) 



n'V 



Here, i^'J'-mJ stands for the probability that a dark atom 
in a state n'V will decay to a state nl. The volumetric 
recombination rate to an atomic state nl is 



OLnliTDM ,Td) = 



(2tt 



,3/2 



(16) 



x / e 
'o 



x [l + fBB(B D (K 2 +n- 2 ),T D )] d{n 2 ). 

Here, n denotes the momentum of the incoming dark 
electron in units of a D /2B D and fBB(E,T) = {e E l T — 
1) _1 is the dark-photon distribution function at energy E 
for a blackbody spectrum at temperature T. The details 
of the atomic physics are encoded in the j n i factor [85] 

7 nJ (/e) = -—^a 3 D B D (l + n 2 K 2 ) 3 



37m 2 



max(l,l')\g(n,l,K,l')\ 2 , (17) 



i'=i±i 



where g(n,l, k,V) denotes the bound-free radial matrix 
elements |92) . Numerically computing the momentum 
integral in Eq. ( 16 ) for each state nl and performing the 
sums in Eq. ( 15 ) up to rt max = 250 yield a recombination 



coefficient of the form 



A d 1 {Tdm, Td ) 



2ira D 



3/2 /rp 

D V 1 DM 



3// 



^2501 £j 



T, 



D 



'-DM ■ 



D 



T, 



D 



(18) 



where G^q are universal dimensionless functions encod- 
ing the details of the atomic physics and its interaction 
with the radiation field. These functions are independent 
of the model parameters (ac, n%D, Bp, £) and therefore 
need to be computed only once. For the purpose of nu- 
merical computation, we tabulate G 2 § and G 2 ,^ on a 
grid of Td/Bd and Tdm /Td values and use an interpo- 
lation scheme to obtain accurate values of the effective 
recombination coefficient^] The photoionization rates 
are related to the above photorecombination coefficients 
through detailed balance 



B D \T D ) 



2ir 



3/2 



DM 



T D ,T D ) 



(19) 



2 The code to compute these coefficients was kindly provided by 
Yacine Ali-Ha'imoud. 



In terms of the universal dimensionless functions, this 
reads 



B D \T L 



*bT D 
3tt 



G 



250 ( u ' !)• 



(20) 



In Fig. [3j we compare the effective total recombination 
rate A d (T D m,T d ) = A^{T DM ,T D ) + A%(T DM ,T D ) 
with the approximate rate given by Eq. ( | 14[ ) . The top 
panel compares the recombination rate when the DM and 
the DR are in thermal equilibrium with Tdm = Td- Wc 
see that the approximate rate ( 14 ) performs reasonably 



well over the temperature range where most of the recom- 
bination is happening (0.007 < T D /B D < 0.02). Most of 
the difference between Eq. ( 14 1 and our exact computa- 



tion can be traced to the fact that the former neglects the 
transitions between the high-n atomic levels and assumes 
that the angular momentum substates are in equilibrium. 
As shown in Ref. [S3], this could in principle be taken 
into account by multiplying Eq. ( 14 ) by a fudge factor 
r~> 1.14. We also show the recombination rate derived 
in Pequignot et al. [51] corrected by this fudge factor as 
used in Recfast [53J. Not surprisingly, this last rate is 
an excellent fit to the exact rate over most of the im- 
portant temperature range. However, we see that both 



the Pequignot et al. and the rate given in Eq. ( 14 1 fails 



at high temperature and to a lesser extent, at low tem- 
perature. For regular atomic hydrogen, these errors are 
inconsequent since most of the baryonic plasma is ionized 
until Tqmb/Bh ~ 0.02 and has mostly recombined be- 
fore Tcmb/Bh ~ 0.003. For weakly-coupled dark atoms 
however, these errors can have a substantial effect on the 
late-time ionization fraction of the dark sector. 

The lower panel of Fig. [3] compares the different re- 
combination rates when the DM temperature differs sig- 
nificantly from that of the DR. This can happen for ex- 
ample for a weakly-coupled dark sector which thermally 
decouples from the DR before the onset of or during re- 
combination. In this case, both the canonical rate ( 14 ) 
and that of Pequignot et al. fail to capture the correct 
temperature dependence over most of the important tem- 
perature range. This should not come as a surprise as 
these two rates are purely functions of Tdm and cannot 
therefore capture the contribution from stimulated re- 
combination which is a function on the DR temperature, 
Td- For regular atomic hydrogen, this effect is unimpor- 
tant since thermal decoupling of baryons happens well 
after recombination. Note that Eq. (14 1 systematically 



underestimates the recombination rate leading to an ar- 
tificially large late-time ionization fraction for weakly- 
coupled dark atoms. Given the sensitivity of the recom- 
bination rate on the DM temperature, it is important 
to accurately capture its evolution through the stage of 
thermal decoupling which we discuss in the next section. 

Without solving any differential equation, it is possible 
to obtain an estimate for the late-time ionized fraction 
of the dark sector by solving the condition 



x D n D {A% + A 2 g)aH, 



(21) 
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here the radiative rate for the Is — > 2p transition to its 
collisional counterpart. The radiative rate is given by 
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The corresponding rate for a dark electron to collisionally 
excite a dark atom from its ground state to the 2p state 
is given by (see e.g. Ref. [95] ) 
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r nM = 0.0177; 
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FIG. 3. Comparison between recombination rates. We have 
chosen the dark sector parameters such that they match those 
of regular atomic hydrogen. We plot the approximate re- 
combination rate given by Eq. ( 14 1 (green short-dashed line) 



as well as our rate computed according to Eq. ( 18 1 includ- 
ing all shells up to n max = 250 (red long-dashed line). For 
comparison, we also show the recombination rate given in 
Ref. |94| corrected by a fudge factor of 1.14 as used in Reef ast 
|93| (black solid line). Top Panel: We compare the rates 
when the DM and DR are in thermal equilibrium such that 
Tdm = Td- Lower Panel: Similar to the top panel but with 
T DM = 0MT D . 



which determines the ionized fraction when the recombi- 
nation process goes out of equilibrium. Using the above 
expression for the recombination rates and evaluating 
them at T D /B D ~ 0.007, we obtain 



x D ~ 2 x 10 



-16 



n D h 2 



-i 



( 



m D 



o.ii J vGcvy Vkevy 



(Bn\ 



(22) 



While this expression is only accurate up to a factor as 
large as 10, it illustrates how the relic ionized fraction 
scales with the dark parameters. 

In our above treatment of dark recombination, we 
have only considered the contribution from radiative pro- 
cesses. In principle, collisional processes could also con- 
tribute to the ionization, recombination, and bound- 
bound rates. To estimate whether we can safely neglect 
the contribution from collisional processes, we compare 



where we have assumed that the dark electrons have a 
Maxwellian velocity distribution with temperature Tdm- 
Demanding that the radiative rate dominates over the 
collisional rate at the onset of dark recombination leads 
to the condition 



m D 
GcV 



n D h 2 

0.11 



> 10 



-10 



(25) 



where we have taken Td ~ Tdm ~ 0.02Bd and xd ~ 1. 
We see that the radiative processes are expected to dom- 
inate over their collisional counterpart unless the dark 
sector is very cold (£ <C 1), dark atoms are very light 
(mD *C 1 GeV) or very weakly coupled. As we will dis- 
cuss in sectionfVl all atomic DM models that are observa- 



tionally relevant obey the condition given by Eq. 25 We 
therefore neglect any contribution from collisional pro- 
cesses in this work. 



C. Thermal Decoupling of Atomic Dark Matter 

In the early Universe, frequent interactions between 
the dark photons and the dark fermions keep the dark 
sector in thermal equilibrium at a single temperature. 
Dark photons Compton scatter off dark electrons, hence 
transferring energy to the DM gas. This energy is then 
redistributed among the dark sector fermions through 
Coulomb scattering between the dark electrons and the 
dark protons. The typical timescale for this process is 



'V-Dl- DM 



\j2-Kct 2 D riDXD In A 
3.8 x KT 4 £5 / 0.11 



(26) 
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where we have assumed Tdm = Td in going from the first 
line to the second equality. This timescale should there- 
fore be considered as an upper limit since after thermal 
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decoupling, we generally have Tdm < Tjj. Here, In A is 
the Coulomb logarithm and is approximately given by 



In A ~ In 



t; 



,3/2 



DM 



3 /2 



(27) 



and has value In A ~ 30 — 60 over the parameter space of 
interest. The Coulomb rate should be compared to the 
Hubble timescale during radiation domination 



rn = \H- 



2.42 x 10 12 
^3.36 + g^ D { 4 



T 

eV 



(28) 



For all interesting atomic DM parameter space, we al- 
ways have T e - P <C Th and therefore, we can always as- 
sume that the dark fermions are in thermal equilibrium 
among themselves at a single temperature T DM . The 
neutral dark atoms maintain thermal contact with the 
dark fermions through elastic collisions. To a good ap- 
proximation, the cross section for a collision between a 
dark electron of energy E e and a dark atom in its ground 
state is given by [9"7] 



320a 2 D 



1 



R2 
D 



(29) 



which is valid for E e < Bd and where 7 ~ 15.69 is 
a best-fit parameter. Averaging this cross-section over 
the Maxwell-Boltzmann velocity distribution of the dark 
electrons and multiplying by the number density of neu- 
tral atoms, we obtain the elastic collision rate between 
dark atoms and dark electrons 



r e 



H 



5.3 x 10 8 a|,(l 
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which is valid for Tdm/Bd < 1. Again, the timescale 
associated with this collisional process r e _H = T^-h is 
much shorter than the Hubble time over most of the dark 
atom parameter space. We can therefore safely assume 
that the whole DM sector (ions + neutral dark atoms) 
is always in thermal equilibrium at a single temperature 
Tdm ■ For the remainder of this section, we thus focus our 
attention on the interaction between the dark fermions 
and the DR. 

Due to its steep dependence on the DR temperature, 
Compton heating is always the dominant energy-transfer 
mechanism between DM and DR at early times. The 
Compton heating rate is given by [72] 



Compton 



me 



64ir 3 a 2 D T D 



xd 



135m 3 , 1 + x£> 



(31) 



where m e and m p are the masses of the dark electron 
and dark proton, respectively. The prefactor in the 



square bracket accounts for Compton heating of both 
dark electrons and protons. For regular atomic hydrogen, 
the large photon-to-baryon ratio ensures that Compton 
heating alone efficiently maintains thermal contact be- 
tween baryons and the photon bath well after the for- 
mer recombine into neutral hydrogen and helium. For 
dark atoms however, there is a large parameter space for 
which Compton heating becomes inefficient at early times 
(i.e. for Tu 3> Bd)- In this case, one must consider other 
possible energy-exchange mechanisms between the DR 
and DM. Before dark recombination, photo-ionization 
heating and photo-recombination cooling are the most 
important mechanisms that can maintain thermal con- 
tact between DM and DR once Compton heating falls out 
of equilibrium. Free-free (Bremsstrahlung) cooling and 
free-free heating are also relevant energy-exchange mech- 
anisms for dark atoms. Finally, dark photons can ex- 
change energy with neutral dark atoms through Rayleigh 
scattering. The volumetric energy- exchange rates for 
these processes are (in energy per unit time per unit vol- 
ume, see Appendix [Cl 



n p _i(7b) 



apT D 
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(34) 



(35) 



where n n 



.Lp-!, -L-i-p-r, Us and are the photo-ionization 
heating, photo-recombination cooling, the net free-free 
heating rate, and the net Rayleigh heating rate, re- 
spectively. Fp—i and F p ^ r are universal dimension- 
less functions parametrizing the details of the atomic 
physics and its interaction with the radiation field. 
These functions are independent of the model parame- 
ters (ap, Bu,mD,i) and therefore need to be computed 
only once. The fractional temperature difference is de- 
noted by e = (T D - T DM )/T D . 

It is instructive to compare the relative magnitude of 
the different energy-exchange mechanisms. The typical 
timescale required by the thermal processes to transfer 
an O(l) fraction of the kinetic energy of the DM is given 
by 



= rr i = ( m * 

\3T DM n D (l + x D ) 



(36) 



where we have assumed the equipartition of energy 
among all dark constituents. Comparing Eqs. (33 1 and 
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FIG. 4. Comparison between the rates of different energy- 
exchange mechanism. We display the rates for Compton heat- 
ing (solid, black), photo-recombination cooling (short-dashed, 
blue), and free- free cooling (dotted, green). We also show the 
Hubble expansion rate (long-dashed, red). The upper panel 
displays the evolution of the thermal rates for an atomic DM 
model with Tbf ~ 700 and T_r ~ 6 x 10 -6 , the middle panel 
has Tbf ~ 280 and T_r ~ 400, while the lower panel shows 
the evolution for a model with Tbf ~ 5 x 10~ 4 and Tr ~ 4. 



(34), we see that the free- free processes and the bound- 
free processes have similar leading-order amplitudes. Ex- 
plicitly taking ratio of the Compton heating rate to that 



of photo-recombination cooling, we obtain 



Compton 



10 



-3 "J 
XD 



Tr 



W\ 5 frno] g 4 (l + -z) 



Br 



(37) 

We first notice that it is always possible to find a high- 
enough redshift such that Compton heating dominates 
over the bound-free energy-exchange channels. As the 
Universe cools down, the photo-heating can become the 
dominant energy-exchange mechanism for low values of 
a.D, for a cold dark sector (£ <C 1), or for a light DM can- 
didate. Generally, photo-recombination cooling, photo- 
ionization heating, and the free-free processes must be 
taken into account when 
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rCompton 
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< 1. (38) 



Here, Zdrcc ^ s ^ ne redshift at which dark atoms recombine. 
For standard atomic hydrogen, we have Tbf ~ 10 2 and 
can therefore neglect any contribution beyond Compton 
heating. 

For relatively large values of the coupling constant a d , 
recombination is generally very efficient, resulting in a 
considerably low ionized fraction xr> < 1 at late times. 
Consequently, thermal coupling mechanisms that depend 
on the presence of free ions to proceed such as Comp- 
ton heating and free-free heating become relatively inef- 
ficient. When this happens, the only remaining mecha- 
nism that can maintain thermal equilibrium between the 
DR and DM for these models is Rayleigh heating. Gen- 
erally, Rayleigh heating takes over Compton heating as 
the dominant heat-exchange channel after dark recombi- 
nation if the condition 



Tr 



rCompton 



Ii 



5 X 10 XD ( itld 
1-Xd 



<1, (39) 



is satisfied. This indicates that only models with very 
small left-over ionization fraction and moderate value of 
the ratio mr>/m e can obtain a non-negligible contribu- 
tion from Rayleigh heating. In Fig. [4j we compare the 
rates for the four dominant energy-exchange mechanisms 
to the Hubble rate. The upper panel shows a relatively 
strongly-coupled dark-atom model with Tr ~ 6 x 10~ 6 
and Tbf ~ 700 where Rayleigh heating becomes the 
dominant thermal-coupling mechanism after the onset of 
dark recombination. In the middle panel, we illustrate a 
model with Tbf ~ 280 and Tr ~ 400 where Compton 
heating is the only important mechanism until adiabatic 
cooling takes over at late times. The lower panel displays 
an alternate scenario with Tbf ~ 5 x 10~ 4 and Tr ~ 4 
where photo-ionization heating and photo-recombination 
cooling dominate from early times until adiabatic cooling 
takes over after recombination. 

We therefore see that the thermal evolution of atomic 
DM strongly depends on the specific choice of dark pa- 
rameters. In contrast to regular atomic hydrogen whose 
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thermal history can be accurately captured by only con- 
sidering Compton heating, it is necessary in general to 
include complementary thermal coupling channels to pre- 
cisely determine the thermal-decoupling temperature of 
dark atoms. Put differently, the thermal history of the 
baryon-photon plasma represents only one possibility 
among all the regimes that a plasma in an expanding 
universe can explore. In this respect, it is interesting to 
realize that the thermal evolution of the baryon-photon 
plasma is rather simple compared to what it could have 
been, had the parameters of SM been different. 

Putting the different pieces together, the Boltzmann 
equation governing the evolution of the DM temperature 
is then l72l 
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(40) 



where the first term of the right-hand side corresponds to 
adiabatic cooling due to the expansion of the Universe, 
the second term takes into account the contribution from 
bound-free and free-free processes as well as Rayleigh 
scattering, and the last term describes the Compton heat- 
ing of the DM gas. 



D. Joint-Evolution of DM Temperature and 
Ionization Fraction 



The thermal history of the dark sector sector is speci- 
fied by simultaneously solving Eqs. Q, and Q 
together with the initial conditions 



XD(zi) = .TD,Saha(^i) 



x 2s (zi) = e 4T o (1 - x D (zi)) 
Tdm{zi) = T D (zi), 



where xn.Saha is the Saha equilibrium ionization fraction 
and z\ is the initial redshift. It is obtained by solving the 
Saha equation 
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These ionization and temperature evolution equations 
are extremely stiff at early times and therefore require 
the use of a stiff solver. We assume that adiabatic cool- 
ing dominates the evolution of the DR temperature such 
that 



T D (z) - T(? MB 



(42) 



where is the temperature of the CMB today. Note 

that dark recombination always happens after the dark 
electrons and positrons annihilate and we have therefore 
neglected the change in the number of relativistic species 
(see Eq. [6| . Equation ( 42 1 is valid as long as we can 



neglect the energy injected into the dark radiation bath 
in a Hubble time, that is 



5p~f,D 
Pf,D 



n„ 



Pj,D 
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(43) 



In the above, we have only included the contribution from 
photo-recombination. In practice, all mechanisms lead- 
ing to a net energy transfer from the DM to the DR 
should be included. Before DM recombination, the right- 
hand side of Eq. ( 43 ) is almost a constant and is equal 
to 



P<y,D 



£ 4 VGeV/ 



Bd 
eV 



< 1. 



(44) 



Models not respecting this bound are likely to require 
a more involved analysis since the DR field cannot be 
taken to be thermal in these cases. A full solution to the 
energy-transfer problem is out of the scope of this paper 
and we therefore focus on models obeying Eq. ( 44 1 . 



We can now compare our improved treatment of the 
dark-atom recombination and thermal decoupling to the 
"standard treatment" described in Ref. [12] • The "stan- 
dard treatment" combines Eqs. ([7]), (JsJ) , and ([£]) into a 
single differential equation for X£> and uses the approx- 
imate recombination rate given in Eq. ( fl4| ). Its DM 
temperature evolution equation only includes Compton 
heating and adiabatic cooling. In Fig. [5j we display the 
ionized-fraction and temperature evolution for a model 
with T R ~ 6 x 10~ 6 and T BF ~ 700. We see that 
the redshift evolution of our improved ionization calcu- 
lation closely matches the standard treatment at early 
times. At late times however, our improved calculation 
predicts an ionization fraction that differs from the stan- 
dard treatment by as much as 60%. The very good agree- 
ment at early times follows from the fact that the ionized 
fraction of the dark sector is initially mostly controlled 
by the Saha equilibrium condition, and therefore insen- 
sitive to the exact value of the recombination coefficient. 
At late times however, the Saha approximation breaks 
down and the ionized fraction becomes sensitive to the 
recombination rates. As can be seen in the lower panel 
of Fig. [5j the inclusion of Rayleigh scattering postpones 
the thermal decoupling of DM, resulting in a DM gas 
that is hotter than one would expect by considering only 
Compton heating. Since the recombination coefficient is 
very sensitive to the DM temperature, this explains the 
somewhat large difference in xu at low redshift. 

Figure [6] displays the ionization and temperature evo- 
lution of a dark-atom model for which Compton heating 
is the dominant thermal coupling mechanism until late 
times (Tbf ~ 2x 10 5 and Tr ~ 10 2 ). In this case, 
the usual calculation accurately captures the behavior of 
both the ionization evolution and the DM temperature. 
This is a not a surprise, since the standard treatment 
was designed to capture this specific regime of the dark 
plasma. The small difference in the ionization fraction 
(up to ~ 12%) at late times is entirely due to our more 
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FIG. 5. Comparison between our improved treatment of 
dark recombination and the standard treatment. We dis- 
play results for a relatively strongly-coupled dark sector with 
Tr ~ 6 x 10 -6 and Tbf ~ 700. The upper panel shows the 
evolution of the ionization fraction as a function of redshift 
while the lower panel shows the corresponding evolution of 
the DM and DR temperatures. 



accurate recombination coefficient which includes the ef- 
fects of high-n shells on the recombination process. 

In Fig. [7J we display the ionization history (upper 
panel) and temperature evolution (lower panel) for a 
model with Tbf ~ 5x 10~ 4 (Tr ~ 4). Not only does 
the standard calculation fail to predict the right recom- 
bination redshift, it also underestimates the late-time 
ionization fraction by more than 50%. Since the Saha 
equilibrium does not hold for these weakly-coupled mod- 
els, their ionization evolution is strongly determined by 
the value of the recombination coefficient. In the lower 
panel, we see that the inclusion of the bound-free and 
free-free processes delays the thermal decoupling of DM, 
hence postponing the onset of dark recombination. As 
the recombination coefficient is very sensitive to the DM 
temperature (see Fig. [3]), this delayed thermal decoupling 
acts to suppress the recombination rate, hence leaving a 
larger ionized fraction at late times. 

In summary, the standard recombination treatment 
originally described in [70] can only be accurately ap- 
plied to dark-atom models for which both bound-free 



FIG. 6. Comparison between our improved treatment of dark 
recombination and the standard treatment. We display re- 
sults for a dark sector where Compton heating always dom- 
inates the thermal coupling between the DR and DM. Here, 
Tbf ~ 1.8 x 10 5 and Tr ~ 10 2 . The upper panel shows the 
evolution of the ionization fraction as a function of redshift 
while the lower panel shows the corresponding evolution of 
the DM and DR temperatures. 



and Raylcigh heating are negligible. For these scenarios, 
the late-time difference between our improved calculation 
and the standard treatment is almost entirely due our 
more accurate recombination coefficient, which properly 
accounts for the effects of excited atomic states and for 
the difference between the DM and DR temperature. On 
the other hand, the standard treatment generally overes- 
timates the thermal-decoupling temperature for atomic 
DM models with Tbf ^ 1 or Tr < 1. As a consequence, 
the standard treatment tends to underestimate the DM 
ionized fraction at late times for these types of models. 



E. Existence of Dark Atoms 

As the dark fine-structure constant is decreased and 
the mass of the dark proton is increased, it becomes 
progressively more difficult for oppositely-charged dark 
fermions to find each other and form neutral bound 
states. There exist critical values of the masses and cou- 
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line corresponds to the bound given in Eq. (451, 
the mostly ionized models (below the line) from the region 
of parameter space dominated by dark atoms at late times 
(above the line). 



FIG. 7. Comparison between our improved treatment of dark 
recombination and the standard treatment. We display re- 
sults for a weakly-coupled dark sector with Tbf ~ 5 x 10~ 4 
and Tr ~ 4. The upper panel shows the evolution of the 
ionization fraction as a function of redshift while the lower 
panel shows the corresponding evolution of the DM and DR 
temperatures. 



pling constants beyond which dark atoms do not form 
and the dark plasma remains ionized even for To <C 
Bjj. Generally, this happens if the recombination rate 
is smaller than the Hubble expansion rate, when it be- 
comes energetically favorable to form dark atoms. Using 
Eqs. (18 1 and (28), dark atoms can form only if the fol- 



lowing condition is satisfied 



n D h 2 
o.n 



/ m D \ - 1 / Bd_ 
VGcV/ UeV 



> 1.5 x 1CT 16 , (45) 



where we have used Tjj ~ Tdm — 0.02Bd, which corre- 
sponds to the usual dark sector temperature at the onset 
of dark recombination. Models violating this bound are 
effectively hidden charged DM models similarly to those 
discussed in Refs. [551 [53] • We display the constraint 
(451 in Fig. [8] (dashed line) together with the values of 
the late-time ionized fraction for a model with Bo = 10 



keV. We see that Eq. (45) delimits very well the region 



where the dark sector is mostly ionized at late times. 



IV. EVOLUTION OF DARK-MATTER 
PERTURBATIONS 

Having determined the evolution of the background 
ionized fraction and temperatures, we now turn our at- 
tention to the evolution of cosmological perturbations in 
the atomic DM scenario. In this theory, the DM and 
DR form a tightly-coupled plasma in the early Universe, 
much like the baryon-photon fluid. Once modes enter the 
horizon, the pressure provided by the relativistic dark 
photons gives rise to a restoring force opposing the grav- 
itational infall of DM, hence leading to dark acoustic os- 
cillations (DAO) in the plasma. Compared to a CDM 
model, the presence of these DAOs delays the onset of 
DM fluctuation growth until the epoch of kinetic decou- 
pling. In addition, atomic DM fluctuations miss out on 
the kick due to the decaying gravitational potential when 
they cross inside the Hubble horizon [45 . We therefore 
generically expect DM fluctuations to be suppressed on 
small scales in atomic DM models. 

We begin this section by giving the key equations de- 
scribing the DM and DR perturbations, emphasizing the 
new collision term between these two constituents of the 
dark sector. We then discuss the different contributions 
to the opacity of the dark plasma and study their im- 
pact on the kinetic decoupling epoch. We then describe 
the various regimes that the perturbations in the dark 
plasma encounter as they evolve through the cosmic ages. 
We finally present numerical examples of these different 
regimes, both in Fourier space and in configuration space. 
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A. Perturbation Equations 

The equations governing the evolution of atomic DM 
fluctuations are very similar to those describing the 
baryon-photon plasma. Special care must however be 
taken to include all the relevant contribution to the opac- 
ity of the DM to DR. The Boltzmann equations for DM 
are 



S D + 9 D - 30 = 0, 



c 2 D k 2 5 D 



k 2 ^ = 



Rd ( 

T~D 



(46) 



%-6 D ), (47) 



where we closely followed the notation of Ref. [HI] in 
conformal Newtonian gauge. Here, Sd is the DM density 
constrast, 8d and 6^ are the divergence of the DM and 
DR velocity, respectively; <j> and ip are the gravitational 
scalar potentials, Rd = Ap^/3pn, cd is the sound speed 
of DM, k is the wavenumber of the mode and r^ 1 is 
the opacity of the dark plasma. Here, the subscript 7 
always refers to the dark photons. The right-hand side 
of Eq. ( 47 ) represents the collision term between the DM 



and the DR. At early times, we generally have Rd 3> 
1 and td <C th, implying that the DM is effectively 
dragged along by the DR. The latter evolves according 
to the following Boltzmann equations: 



+ -0y - 4(j> = 0; 



(48) 
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Eqs. (48) and (49) describe the evolution of the dark- 



photon over-densities (&y) and of the dark-photon veloc- 
ity, respectively. It is also necessary to solve for the hier- 
archy of dark-photon multipoles (Eqs. (50) and (51)) to 
properly account for DR diffusion and its impact on DM 
perturbations. Since our focus is to describe the clus- 
tering of DM in this model, we do not solve for the DR 
polarization. 

During the radiation-dominated epoch, the energy 
density of the DR is generally subdominant compared to 
the contribution from regular photons and neutrinos (see 
Eq. ([5])). Therefore, the time-dependence of the gravita- 
tional potentials <f> and ip is very similar to the standard 
CDM case for which we have 



-3c/>p 



sin (kr/VS) - (kr/y/3) cos {kr/y/S) 
(fcr/v^) 3 



, (52) 



where r stands for the conformal time and cp p is the pri- 
mordial amplitude. Since the gravitational potential is 
an oscillatory function, Eqs. (|46|) and ([48]) essentially de- 



scribe driven harmonic oscillators where the driving force 
is provided by the baryon-photon plasma. Indeed, taken 
as a whole, the equations describing the dark plasma and 
the baryon-photon plasma in the early Universe corre- 
spond to a system of coupled harmonic oscillators. Such 
a system is known to exhibit resonance phenomena when- 
ever the driving frequency approaches the natural fre- 
quency of the oscillator. For the dark plasma however, 
we always have 



R D « R = 



3pb 



(53) 



where p 7 and pb are the energy densities of regular pho- 
tons and baryons, respectively. As the sound speed of a 
plasma is approximately given by c p — 1/ \/3(l + i? _1 ), 
Eq. ( 53 ) implies that the sound speed of the dark plasma 



is always smaller than that of the baryon-photon plasma. 
Thus, the dark plasma is never driven close to its reso- 
nance threshold^] 



B. Dark Opacity and Kinetic Decoupling 

The opacity of the dark plasma dictates the strength of 
the coupling between DM and DR. Heuristically, td can 
be considered as the mean free path a dark photon trav- 
els between collisions with a dark ion or a dark atom. As 
such, Rd/td is approximately the momentum-transfer 
rate between the DR and the DM. As in the case of the 
thermal coupling of DM to DR, many mechanisms con- 
tribute to the exchange of momentum between the two 
dark components. In addition to the usually-considered 
Compton-scattering term, we also include the contribu- 
tion from Rayleigh scattering as well as the contribution 
from bound-free processes. The opacity is then 



' D ~ ' Compton ' R T 'p-i 

The Compton scattering term is given by 



Compton 



an D XD<JT,D 



1 



(54) 



(55) 



where ctt,_d = 8Tra D /(3m 2 ) is the dark Thomson cross 
section and a is the scale factor. The factor in the bracket 
accounts for Compton scattering off dark protons. The 
contribution from Rayleigh scattering can be written as 



r R 1 = an D (l - Xd)(<Jr) 

~ 327T 4 an£>(l — xd)<tt,d 



Td_ 
B d 



(56) 



3 While resonant enhancement does not occur for the simple 
atomic DM scenario considered in this paper, it is nevertheless 
possible to construct a model where such resonance happens. 
This is an interesting possibility that we leave for future work 
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which is valid for To <C Bp . Here ctr is the cross-section 
for Rayleigh scattering (see Eq. (C13)) and (. . .) denotes 



a thermal average with respect to the Planck function 
describing the distribution of dark photons. Finally, the 
photoionization contribution is 



Tp- 1 ! = a n m (&nl 



(57) 



an D x 2s 



6C(3)T£ 



D_ -B D /(4T D ) a 2s ( T^D_ i 



5d 



where n n i and cr„; are the number density and the pho- 
toionization cross-section for dark atoms in quantum 
state nl. 

For most atomic DM models, the Compton scatter- 
ing contribution to the opacity dominates before dark 
recombination, with photoionization giving a subdomi- 
nant contribution for models with Tbf ^ 1- Once a sig- 
nificant number of dark atoms has recombined, Rayleigh 
scattering can become the dominant source of opacity for 
models with Tr < 1. 

The DR effectively begins free-streaming when its 
mean-free path becomes comparable to the size of the 
Hubble horizon, that is, r^ 1 ~ H. In a model for which 
Compton scattering dominates the interactions between 
the DR and the DM, the onset of the free-streaming 
epoch is independent of the temperature and only de- 
pends of the fraction of ionized dark atoms, 
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(58) 
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where we have assumed that the DR decoupling happens 
in the radiation-dominated epoch and have neglected the 
contribution from dark protons. If the second argument 
in the min [.., ..] function is larger than unity, it indicates 
that dark photons begin free-streaming at a very early 
epoch, making them difficult to distinguish from massless 
neutrinos. For a strongly-coupled dark sector which have 
^5 1) the decoupling of dark photons happens much 
later and is governed by Rayleigh scattering of dark pho- 
tons off neutral dark atoms. This leads to a decoupling 
temperature 
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7 x 10" 
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(59) 

where we have taken the limit xu|dec « 1 It is impor- 
tant to note that for most atomic DM scenarios, Comp- 
ton scattering alone is responsible for the DM-DR cou- 
pling. 

The key quantity governing the size of the smallest DM 
structure in the Universe is the temperature at which 
DM kinematically decouples from the DR (i.e. the dark 
drag epoch). This temperature can approximately be 



determined by solving the condition Rd/td ~ H. For 
models dominated by Compton scattering, this leads to 
a kinetic decoupling temperature given by 
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(60) 



This is effectively an implicit equation for the drag-epoch 
temperature since the right-hand side involves the ionized 
fraction of dark atoms evaluated at that epoch. Taking 
££>|drag ~ 1 leads to a lower bound on the temperature 
at which DM ceases to be dragged along by the DR. If 
< 1) the kinetic decoupling occurs when Rayleigh 
scattering becomes ineffective. In this case, we obtain 
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B D 
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(61) 



where we have used x £>|drag *C 1, which is a necessary 
condition for Rayleigh scattering to dominate over Comp- 
ton scattering. 

DM fluctuations on subhorizon length scales at the 
drag epoch have a significantly different evolution than 
those of a standard ACDM model. In particular, the 
suppression of small-scale power due to the dark-photon 
pressure leads to a minimal halo mass at late times. We 
will discuss this effect further in section IV Al 



C. Regimes of the Dark Plasma 



We now turn our attention to the formal solutions to 



Eqs. (461 to (51 1. These admits different regimes depend- 



ing on the relative values of the opacity, wavenumber, and 
Hubble expansion rate. These regimes are 

1. Superhorizon Regime: This regime occurs when 
the wavelength characterizing a mode is still larger 
than the Hubble horizon, k < H. As for the 
CDM case, cosmological perturbations do not sig- 
nificantly evolve in this regime. We therefore do 
not discuss this case any further since atomic DM 
is indistinguishable from standard CDM for these 
modes. 

2. Dark-Acoustic-Oscillation (DAO) Regime: 

Once modes cross into the Hubble horizon, the mi- 
crophysics governing the interaction between the 
DM and the DR becomes effective. If the mean- 
free-path of dark photons between collisions with 
dark fermions is much smaller than the wavelength 
of a given mode (kru C 1), then we can consider 
the DM and the DR to be tightly-coupled and to 
form an almost perfect fluid at this scale. In this 
case, the DR pressure can effectively counteract 
the pull from the gravitational potential, leading 
to the propagation of acoustic oscillations in the 
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dark plasma. This DAO regime occurs for wave- 
lengths larger than the diffusion length scale of the 
dark photons, that is, 



H < k < k 



D 



(62) 



where kr> is diffusion damping scale and is defined 
by [98] 



1 



k 2 D (r) 



T D dr' 



6(1 + R} 



1 



Rji + Rd 



(63) 



Heuristically, 1/ko corresponds approximately to 
the average distance travelled by a dark photon in 
a Hubble time. Note that Eq. ( 63 ) is the result of a 



first-order expansion in krr> of the dispersion rela- 
tion of DAOs and therefore is only accurate in the 
limit krn C 1. In the DAO regime, DM fluctua- 
tions undergo constant-amplitude oscillations 



Sd ~ exp <ik drc 



(64) 



where c p = 1/\J 3(1 + R^) is the sound speed of 
the dark plasma. Due to these oscillations, DM 
fluctuations entering the horizon in this regime miss 
out on both the logarithmic growth and the horizon 
kick. Thus, despite the absence of damping in the 
DAO regime (besides the small contribution from 
Hubble expansion), DM perturbations still gener- 
ally display smaller amplitudes at these scales when 
compared to a standard CDM model. 

Diffusion-Damping Regime: Once the average 
distance travelled by a dark photon during a Hub- 
ble time becomes comparable with the wavelength 
of a given mode, DR can effectively diffuse out of 
over-densities at these scales. If the momentum 
transfer rate between the DR and the DM is larger 
than the plasma oscillation frequency 



k D < k < 



±Rd 

Cp To ' 



(65) 



then DM is dragged along by the escaping dark 
photons, effectively erasing DM fluctuations on 
these scales. The time-dependence of the DM per- 
turbations is then given by 



8d ~ exp lik I drdp > exp 



w 



(66) 



leading to an exponential damping on fluctuations 
on scale smaller than the diffusion length. As be- 
fore, we have assumed the dark plasma is tightly- 
coupled (krn < 1) in deriving Eq. (66). When 



this condition ceases to be satisfied, the damping 
of DM fluctuations is no longer exponential and 
Eq. (66) breaks down. Note that the diffusion 
damping regime ends when kc p ~ Rd/td since DM 
ceases to be dragged by the DR when this condition 
is satisfied. 



4. Acoustic-Damping Regime When the mean- 
free-path of dark photons becomes comparable to 
the wavelength of a given mode (fcr^ > 1), the 
dark plasma ceases to behave like a single perfect 
fluid. Due to the slow reaction of DM to the mo- 
tion of dark photons in this regime, acoustic os- 
cillations in the DM fluid start to significantly lag 
those propagating in the DR fluid. This in turns 
leads to viscous dissipation in the dark plasma re- 
sulting in the damping of DM fluctuations at these 
scales. Note however that this acoustic damping is 
weaker than the exponential damping characteristic 
of the tightly-coupled diffusion regime. Generally, 
acoustic damping occurs for modes satisfying the 
condition 



TD C p T D 



(67) 



The last inequality ensures that a significant 
amount of momentum is still transfered to the DM 
during acoustic damping. Physically, this regime 
is characterized by the development of anisotropic 
stress in the dark fluid associated with a significant 
dark-photon quadrupole moment. As such, it is 
therefore difficult to obtain an analytical solution 
to the evolution of DM fluctuations in this regime. 

Gravitationally-Dominated Regime: When 
the momentum transfer rate between the DR and 
atomic DM falls below the oscillation frequency of 
the dark photons, that is 



H < — < kc„ 

TD 



(68) 



the dark plasma ceases to behave like a single fluid 
and we must consider the DM and the DR fluctu- 
ations separately. In this regime, the evolution of 
the DM perturbations is determined by a competi- 
tion between the gravitational potential dominated 
by the baryon-photon plasma and the remaining 
pressure of the DR. Since both of these contribu- 
tions are oscillatory in nature, the evolution of the 
DM fluctuations in this regime can be rather com- 
plex. Essentially, we can view DM perturbations 
as forced oscillators driven by two distinct forces 
oscillating with difference frequencies. Explicitly, 
the equation governing the evolution of DM fluctu- 
ations is 



<h. ( # + — W = S G (k,r) - (69) 

TD J T D 



where the gravitational source term is given by 
S G (k,T) = 3$ - *V + 3^ + (70) 
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The dark-photon driving term can be obtained 
from rearranging Eqs. (49) and (48) 



7ty — od- 



3t d 
~ 4 



k 2 k 2 



(71) -3 



In the limit that a DM fluctuation enters this 
regime after being exponentially damped by the 
diffusion of dark photons, the evolution of DR per- 
turbations is dominated by the gravitational po- 
tential. Neglecting the small quadrupole moment 
of the dark photons, this implies that the term in 



the square bracket of Eq. (|71| is very close to z ero 
We thus have 6^ 
reduces to 



0D in this limit and Eq. (69) 



S D + HS D = 3<f> - k 2 i/j + ZH<t>. 



(72) 



Eq. ( 72 ) is to the usual equation describing the 
growth of standard CDM fluctuations. It implies 
that DM fluctuations grow logarithmically with the 
scale factor during the radiation-dominated epoch. 

On the other hand, for modes that enter the 
gravitationally-dominated regime before being sub- 
stantially damped by photon diffusion, the dark- 
photon driving term cannot be neglected in 
Eq. ( 69 ) . In this case, the dark-photon quadrupole 



moment is non- negligible in Eq. (71), leading to an 
oscillatory pressure term with 



e ^ 



(73) 



The gravitational potential contribution Sc(k,T) 
also oscillates, but with a somewhat lower fre- 
k 



quency ~ k/y/3(l + R). As these oscillators 
can randomly go in-phase and out-of phase, some 
damping or even some amplification can occur in 
this regime. However, since the momentum transfer 
rate Rd/td is decaying more rapidly than Sa(k,T), 
the gravitational potential rapidly becomes the 
dominant contribution to the evolution of DM fluc- 
tuations. Indeed, while Rd/td oc t~™ with n > 3 
(depending on which process dominates the calcu- 
lation of the opacity), Sc{k,T) oc r~ 2 . Therefore, 
the DM fluctuations rapidly settle into the loga- 
rithmic growing mode. 

Kinetically-Decoupled Regime: As the mo- 
mentum transfer rate between atomic DM and DR 
falls below the Hubble rate, Rd/td < H, DM es- 
sentially stops interacting with the dark photons 
and begin behaving like the standard CDM. Modes 
entering this regime during radiation domination 
begin growing logarithmically with the scale fac- 
tor while modes entering during matter domination 
grow linearly with a. In all cases, modes that en- 
ter the Hubble horizon after kinetic decoupling are 
undistinguishable from those of standard CDM. 
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FIG. 9. Redshift evolution of an atomic DM fluctuation with 
k — 70 Mpc -1 in synchronous gauge (green solid line). We 
identify on the figure the different regimes that the DM fluc- 
tuation encounters during its evolution. These are: (a) super- 
horizon; (b) DAO; (c) diffusion damping; (d) gravitationally- 
dominated; and (e) kinetically-decoupled. For the case dis- 
played here, the fluctuation enters the horizon while the dark 
plasma is tightly-coupled and therefore undergoes acoustic 
oscillations. Once ko ~ k, the fluctuation becomes exponen- 
tially damped by dark-photon diffusion. This damping ceases 
to be effective when the rate of momentum transfer between 
the DR and DM falls below the oscillation frequency of dark- 
photon fluctuations. When this happens, the fluctuation en- 
ters the gravitationally-dominated regime where the pertur- 
bation rapidly settles in tothe logarithmic growing mode. The 
sharp feature around (z + l) -1 = 3 x 1CP 5 is an artifact of 
plotting the absolute value of the fluctuation. For compari- 
son, we also show the behavior of a CDM fluctuation with the 
same wavenumber. 



D. Numerical Solutions 



We solve numerically Eqs. ( 46|[5l | together with the 
standard Boltzmann equations describing the evolution 
of baryons, photon and neutrinos [5T]. We use a modified 
version of the publicly-available code CAMB [99] assuming 
that all of the DM is made of dark atoms. We modify 
the perturbed Einstein equation to include the new con- 
tributions from DM and DR. We use a flat background 
cosmology compatible with WMAP 7- year release jlOOJ : 
il b h 2 = 0.0226, n D h 2 = 0.1123, H = 70.4km/s/Mpc, 
A|, = 2.3 x 10~ 9 , n s = 0.963, and r ro = 0.088. We 
consider pure adiabatic initial conditions 



S D (zi) = 6 b (zi) 6~(zi) = 6Jzi), 



9 D {zi) = Oj(zi) = 9 1 {z- 1 ), 



F ¥ =0, I > 2. 



(74) 



(75) 



(76) 



where z\ is the initial redshift which is determined such 
that all modes are superhorizon at early times, hr(z\) <C 
1. We first pre-compute the ionization and thermal his- 



tory of the dark sector as described in section III and use 
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FIG. 10. Redshift evolution of an atomic DM fluctuation with 
k = 70 Mpc -1 in synchronous gauge (green solid line). We 
identify on the figure the different regimes that the DM fluc- 
tuation encounters during its evolution. These are: (a) super- 
horizon; (b) DAO; (c) Acoustic damping; (d) gravitationally- 
dominated; and (e) kinetically-decoupled. For the case dis- 
played here, the fluctuation enters the horizon while the dark 
plasma is weakly-coupled and therefore transitions to the 
acoustic damping regime very rapidly. This damping ceases 
to be effective when the rate of momentum transfer between 
the DR and DM falls below the oscillation frequency of dark- 
photon fluctuations. When this happens, the fluctuation en- 
ters the gravitationally-dominated regime where the pertur- 
bation is further damped before it settles into the logarithmic 
growing mode. For comparison, we also show the behavior of 
a CDM fluctuation with the same wavenumber. 



the result to compute the opacity of the dark plasma as 
given in Eq. (54) above. The ionization history of the 



baryon-photon sector is pre-computed in the usual way 
using RecFast [93 . The linear perturbation equations 
are then evolved forward in time from z\ to z = 0. At 
early times when krr> <C 1 and td/t <C 1, Eqs. (47) 



and ( 49 1 are very stiff and we use a second-order tight- 



coupling scheme similar to that used for the baryon- 
photon plasma at early times [101H103] . 

In Fig.|9j we show the time evolution of a single Fourier 
mode for relatively strongly-coupled dark atoms. We 
clearly identify the different regimes that the fluctua- 
tion encounters from its horizon crossing to its late-time 
growth. For the particular choice of parameters displayed 
here, this Fourier mode enters the horizon in the DAO 
regime and oscillates until dark-photon diffusion expo- 
nentially suppresses its amplitude. Once DR kinemati- 
cally decouples from DM, the fluctuation can start grow- 
ing like regular CDM but from a much-reduced ampli- 
tude. 

To contrast, we show in Fig. [10] the redshift evolution 
of the same Fourier mode for a weakly-coupled model 
of dark atoms. In this case, kro ~ 1 shortly after the 
mode enters the horizon and therefore it only briefly ex- 
periences the DAO regime. It then rapidly transitions 
to the acoustic damping regime where its amplitude de- 
cays, though not as quickly or strongly as in the diffusion 



damping regime. After kinetic decoupling, the mode fi- 
nally settles into the logarithmic growing mode. 



E. Real-Space Evolution: DAO Scale 

The Fourier-space description of cosmological fluctua- 
tions allows one to qualitatively understand the different 
stages of their evolution and to obtain accurate numer- 
ical solutions to their equations of motion. Ultimately 
however, physical density fluctuations evolve in configu- 
ration space and it is thus important to translate their 
Fourier-space behavior into this latter space. 

In Fig. [Tl] we display the configuration-space redshift 
evolution of a Gaussian adiabatic density fluctuation. 
In panel (a), we see both a dark-plasma and a baryon- 
photon plasma density wave traveling outward from the 
initial over-density. Since the dark plasma generically 
has a lower sound speed than the regular baryon-photon 
plasma, density waves propagating in the former do not 
travel as far in a given time interval as waves propagat- 
ing in the latter. This results in the dark-plasma den- 
sity wave lagging behind its baryon-photon counterpart. 
Panel (a) of Fig. 11 also clearly shows the damping of the 
initial density fluctuation at short length-scales resulting 
from the outward propagation of acoustic waves. 

Panel (b) presents a snapshot of the outward-moving 
waves shortly after dark photons cease to be tightly- 
coupled to the DM. Panel (c) further shows the DR dif- 
fusing out of the DM fluctuation after they kinematically 
decouple from each other. After the dark decoupling 
epoch, DM fluctuations begin to grow while DR contin- 
ues to propagate out of the initial overdensity, eventually 
overtaking the baryon-photon sound horizon. Note that 
the sound horizon of the dark plasma remains imprinted 
on the DM fluctuations, resulting in a preferred scales in 
the late-time density field similar to the standard BAO 
scale. This so-called DAO scale constitutes a tell-tale sig- 
nature of the presence of a dark plasma distinct from the 
standard baryon-photon plasma in the early Universe. 
We will discuss its cosmological implications in the next 
section. 

The inset inside panel (d) shows that the dark-photon 
density wave can actually overtake the baryon-photon 
plasma wave after the former decouples from the DM. 
As in the case of free-streaming neutrinos, this effectively 
establishes gravitational-potential perturbations beyond 
the sound horizon of the baryon-photon plasma. It 
has been shown |104) that such supersonic gravitational- 
potential perturbations can have a measurable impact on 
the CMB power spectrum. This indicates that DR fluc- 
tuations could in principle have an impact on the CMB, 
which we consider section fVl 



Panels (e) and (f) of Fig. 11 show the late-time decou- 
pling of photons and baryons. As in the standard CDM 
scenario, baryons fall in the gravitational potential wells 
established by DM once they cease to be dragged along 
by the photons. In the atomic DM scenario however, the 
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FIG. 11. Redshift evolution of an initially Gaussian-shaped density fluctuation in configuration space. We take the fluctuation 
to be adiabatic. Here, an = 0.08, Bd = 1 keV, mo = 1 GeV, £ = 0.37, and r stands for the comoving spatial separation. 
We display the evolution of the DM (dot-dashed, blue), DR (dashed, green), baryons (dotted, red), and photons (solid, black). 
Panel (a) shows the initial Gaussian density fluctuation together with four snapshots of the outgoing density waves at lower 
redshifts. The next five panels display different stages of the evolution of the dark and baryon-photon plasmas. Note the 
changing axes from panel to panel. The insets in panels (d), (e) and (f) focus on the progression of the baryon-photon sound 
horizon. At late times, two key lengthscales emerge: the standard BAO scale at tbao — 147 Mpc and the new DAO scale at 
tdao — 3.3 Mpc (for the model plotted here). 
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DAO scale becomes imprinted on the baryons as they fall 
toward DM in a process similar to how the BAO scale is 
imprinted on the DM as it falls toward the baryons at 
late times. 

In summary, compared to a vanilla ACDM cosmology, 
the atomic DM scenario has an additional cosmological 
length scale corresponding to the size of the dark-plasma 
sound horizon when the DR kinematically decouples from 
the DM. Since the dark plasma decouples earlier and has 
a lower sound speed than the baryon-photon plasma, this 
DAO scale is generically much smaller than the BAO 
scale. Furthermore, the amplitude of density fluctuations 
on length scales smaller than the DAO scale is severely 
suppressed, leading to a minimal mass for the first objects 
that collapse at late times. We discuss the cosmological 
implications of atomic DM in the next section. 



V. COSMOLOGICAL IMPLICATIONS 

The atomic DM scenario alters cosmological observ- 
ables in two possible ways. On the one hand, the pres- 
ence of dark relativistic degrees of freedom modifies the 
cosmological expansion history of the Universe. On the 
other hand, the gravity and pressure of the dark-photon 
perturbations impacts the evolution of DM and baryon- 
photon fluctuations, hence affecting late-time observables 
such as the CMB and the matter power spectrum. Unfor- 
tunately, the former case does not lead to unique effects 
on observable cosmological probes since the background 
cosmology of atomic DM is indistinguishable from that 
of a standard ACDM model that contains extra relativis- 
tic species. As the impact of relativistic neutrinos on 
the CMB and the matter power spectrum has been ex- 
tensively studied in the literature (see e.g. |104[ 1105] ). 
we shall only briefly review their key effects in the two 
following subsections. 

Conversely, dark-photon perturbations can affect cos- 
mological observables in a way that distinguish them 
from ACDM models containing extra neutrinos. This 
stems from the fact that relativistic neutrinos can free- 
stream from a very early epoch while dark photons can 
only do so after they decouple from the DM. Cosmolog- 
ical modes entering the horizon while the dark photons 
are free-streaming are expected to behave similarly to a 
ACDM model with an equivalent number of relativistic 
neutrinos. On the other hand, modes crossing the hori- 
zon while the dark photons are tightly-coupled to the DM 
do not experience the damping and phase shift of acous- 
tic oscillations 104 usually associated with the presence 
of extra radiation. Therefore, we expect the atomic DM 
scenario to leave a distinct imprint on the CMB if dark 
photons begin free-streaming when the length scales rel- 
evant for this cosmological probe are crossing into the 
Hubble horizon. A corollary of this statement is that a 
dark-atom model for which dark photons decouple very 
early is, as far as the CMB is concerned, indistinguishable 
from a ACDM universe containing extra neutrinos. 



Atomic DM itself alters cosmological observations 
through the modified growth of its density fluctuations. 
As explained in section [TV C[ DM fluctuations entering 
the horizon before kinetic decoupling in the radiation- 
dominated era miss out on the amplitude boost due 
to the rapidly decaying gravitational potential. Fur- 
ther, atomic DM perturbations on scales smaller than 
the dark-photon diffusion length are strongly damped, 
leading to an absence of cosmological structure at these 
scales. The acoustic oscillations in the DM plasma which 
are one of the key feature of the atomic DM scenario 
remain imprinted on the late-time matter power spec- 
trum. In configuration space, these oscillations point 
to an important length scale, the DAO scale, at which 
the clustering of DM is enhanced and below which it is 
suppressed when compared with an equivalent ACDM 
model. For most atomic DM models that are in agree- 
ment with current observations, the DAO scale must lie 
today in the highly non- linear regime of cosmological fluc- 
tuations, hence making its impact on cosmology difficult 
to observe. 



A. Matter Power Spectrum 

The most dramatic and distinct cosmological impli- 
cation of the atomic DM scenario is the modification 
of the small-scale matter power spectrum. Ultimately, 
this is a consequence of the relatively late kinetic decou- 
pling of atomic DM compared to a standard WIMP CDM 
model. Indeed, the delayed kinetic decoupling of the dark 
plasma considerably impedes the growth of DM fluctu- 
ations for all subhorizon modes. Furthermore, fluctua- 
tions on length scales shorter than the diffusion distance 
of dark photons are exponentially damped compared to 
a CDM model, effectively prohibiting the formation of 
any late-time structure at these scales. Subhorizon scales 
that exceed the diffusion damping scale display acoustic 
oscillations that remain clearly imprinted on the late- 
time matter power spectrum, since DM forms the bulk of 
the non-relativistic matter. Fluctuations on length scales 
that cross into the horizon after DM kinematically de- 
couples have a growth history similar to that of a vanilla 
ACDM model and we therefore expect no particular sig- 
nature at these scales (see however a caveat at the end 
of this subsection). 

In Fig. |12| we show examples of late-time linear matter 
power spectra for different models of atomic DM. The 
features discussed above are clearly visible in the spectra. 
We also observe that as the binding energy of the dark 
atoms is increased, the comoving wavenumber at which 
the power spectrum significantly departs from the ACDM 
case is increased. This makes sense since a higher binding 
energy implies an earlier dark recombination and kinetic 
decoupling, therefore pushing the impact of the DM-DR 
coupling toward higher comoving wavenumbers. 

In the two lowest panels of Fig. [l2j we keep the dark- 
atom binding energy and mass constant but vary the dark 
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FIG. 12. Total linear matter power spectrum at z — for 
three atomic DM models. For reference, we also display the 
linear matter power spectrum for a vanilla ACDM model. 



fine-structure constant ap. We observe that as ajj is 
decreased, the smallest comoving wavenumber affected 
by the DAOs moves toward higher values. At first, this 
seems counterintuitive as a higher value of the dark cou- 
pling constant generally leads to a lower residual ioniza- 
tion fraction which in turns allows the dark photons to 
rapidly decouple from the DM. We must however remem- 
ber that dark photons also interact with neutral dark 
atoms through Rayleigh scattering. For a fixed binding 
energy and DM mass, the Compton-scattering contribu- 



tion to the dark-plasma opacity is roughly independent 
of ao since xd{Td <C -Bd) oc a^ 6 and ctt,d a %- 
On the other hand, the Rayleigh scattering contribution 
to rrj is a steep function of old with 1 oc a D after 
the onset of dark recombination. Thus, an increase of 
the dark fine-structure constant can considerably boost 
the Rayleigh-scattering contribution to the opacity of the 
dark plasma, hence significantly postponing its epoch of 
kinetic decoupling. 

It is also instructive to consider the correlation func- 
tion of matter fluctuations in configuration space. The 
linear correlation function is related to the matter power 
spectrum via a 3D Fourier transform which, after simpli- 
fication, can be reduced to 



a(r) 



1 

27r2 



dk k 2 P(k)j (kr), 



(77) 



where Jo (At) is the Spherical Bessel function of order 
0. In Fig. |13| we display the linear correlation function 
computed from the three matter power spectra shown in 
Fig. 12 as well as the correlation expected from a stan- 
dard ACDM model. In all cases, the usual BAO scale 
is clearly visible around r ~ 104/i Mpc. In a similar 
manner, the novel DAO length scale appears as a local 
enhancement of the correlation function at the scale cor- 
responding to the sound horizon of the dark plasma at 
dark decoupling. On scales smaller than this sound hori- 
zon, the correlation function is significantly damped com- 
pared to the ACDM case, a consequence of the damping 
of small-scale fluctuations discussed above. 

At late times, the key signature of these new features 
in the matter power spectrum and correlation function is 
a minimal DM halo mass. Indeed, since most of primor- 
dial fluctuations on scales smaller than the dark-plasma 
sound horizon are effectively wiped out by the diffusion 
of dark photons, no self-bound object can form at late 
times at these scales. The first regions that can collapse 
into self-bound DM halos must then have a minimal ini- 
tial comoving size ~ tdao- Therefore, the first DM halos 
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FIG. 13. Linear correlation function for the three atomic DM 
models plotted in Fig. [12] 
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have a minimal mass given approximately by 



4?r 



r DAO^£ ) Pcrit- 



(78) 



In the hierarchical model of structure formation, these 
first halos are then accreted into larger, more massive 
halos. While some of these minimal-mass halos are de- 
stroyed through tidal stripping in larger halos, it is ex- 
pected that a certain fraction of them remains as discrete 
sub-halos in larger bound objects such as galactic ha- 
los. These could potentially be detected through strong- 
lensing studies of substructures inside galactic halos (see 
e.g. [106 ). Indeed, there is currently a growing scien- 
tific effort aimed at developing analytical methods and 
experimental techniques geared toward the detection of 
small-scale substructures inside galactic halos |107H113] , 
While current results mostly point out the existence of in- 
dividual massive subhalos in galactic strong lenses [1141 - 
I116j . statistical analyses of multiple- image lenses have 
the potential to lead to strong constraints on the possi- 
ble minimal subhalo mass. Thus, it is not unreasonable 
to think that this tell-tale signature of atomic DM might 
be detected in the near future. 

Currently, the most stringent constraints on the small- 
scale matter power spectrum come from the Ly-a for- 
est data [117H119] . There is however considerable sys- 
tematic uncertainties in converting the Ly-a flux power 
spectrum to the actual linear matter power spectrum 
j!20j . Most studies have assumed a power-law spectrum 
with a running spectral index. Since atomic DM pre- 
dicts a much more complex shape of the matter power 
spectrum, it is not straightforward to apply these con- 
straints to this scenario. In reality, hydrodynamical sim- 
ulations of the Ly-a flux power spectrum in the atomic 
DM scenario will be required to derive the appropriate 
Ly-a constraints. We can nevertheless use the current 
measurements to obtain rough guidelines. Ref. [121] has 
found no deviation from the standard ACDM scenario 
on scales k < 2 h Mpc" 1 , while Refs. PT3 Oil Ell] have 
determined that the linear power spectrum is consistent 
with CDM on scales k < 5ft Mpc -1 . Pending unforeseen 
physical effects, it is unlikely that atomic DM can have a 
large impact on these scales while leaving the Ly-a flux 
power spectrum unchanged. We therefore demand that 
the atomic DM linear power spectrum does not signifi- 
cantly deviate from that of a vanilla ACDM cosmology on 
those scales. More precisely, we compute two constraints: 
koAO > 1ft Mpc -1 and /cdao > 5ft Mpc -1 , where we 
have defined fcrjAO = 7r / r DAO- These bounds correspond 
to M min fa 10 13 M Q and M min pa 9.3 x 10 10 M Q , respec- 
tively. 

We show this constraint in Fig. [20] where we display 
countours of constant minimal halo mass (dotted white 
lines) in the ac-mj) plane for various values of the atomic 
binding energy. The dark and light gray regions corre- 
spond to M min > 10 13 M Q and M min > 9.3 x 1O 1O M , 
respectively. We see that the Ly-a forest data only con- 
strain models with Bo < 10 keV. For higher values of 
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FIG. 14. Dimensionless linear matter power spectrum at z = 
for a single dark-atom model for different values of £. We 
plot k 3 P(k) to magnify the small-scale region of the spectrum. 
Throughout the plot, we fix z eq , Q.bh 2 , and the angular scale 
of the CMB sound horizon at decoupling (6 S ). 



the atomic binding energy, the kinetic decoupling of the 
DM happens very early (even for ~ 0(1)) and leads 
to a minimal halo mass M min < 10 7 M which is uncon- 
strained by data. 

To modify the faint-end of the galaxy luminosity func- 
tion and bring it in agreement with the data, atomic 
DM needs to modify the properties of halos in the range 
1O 8 M < M min < 1O 1O M gD]. According to the 
constraints in Fig. 20 this is possible for model with 
100 TeV and a D ~ 0.2. Detailed 



B D ~ 10 keV, m D 
Af-body simulations will be required to assess how the 
atomic DM scenario exactly affects the halo mass func- 
tion and density profile, but it is clear that there are 
allowed values of the dark parameters that can directly 
address some of the dwarf-galaxy problems. 

For completeness, we also mention in passing that the 
dark radiation has a small impact on the evolution of DM 
perturbations even for modes that enter the horizon af- 
ter dark kinetic decoupling in the radiation era. Indeed, 
the free-streaming dark radiation amplifies the DM fluc- 
tuations as they enter the horizon due to their impact 
on the rapidly decaying gravitational potential. To illus- 
trate this effect, we must however be careful since the 
presence of extra radiation shifts the epoch of matter- 
radiation equality z cq and the angular scale of the CMB 
sound horizon at decoupling 9 S . Since these quantities 
have been precisely measured by experiments, we should 
keep them fixed as we vary the amount of dark radiation. 
This entails to adjusting the physical DM density (f^/i 2 ) 
and the Hubble constant to keep the redshift of equality 
and the angular scale at decoupling unchanged. We keep 
ilfeft 2 fixed throughout since the CMB tightly constrains 
its value. 



We illustrate in Fig. 14 the linear matter power spec- 
trum for different values of the ratio of the dark sector 
temperature to the CMB temperature. We clearly ob- 
serve that as £ is increased, the amplitude of Fourier 
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modes crossing into the Hubble horizon in the radiation 
era are enhanced. In addition to the previously men- 
tioned effect caused by the decaying gravitational poten- 
tial, this enhancement is also caused by the larger DM 
density of model with high values of £. Indeed, the mat- 
ter fluctuation spectrum is sensitive to the ratio n&/f2 m , 
since baryons are withheld from gravitational collapse 
prior to hydrogen recombination. As the total matter 
density is increased for a fixed amount of baryons, more 
matter can form gravitationally-bound structure and the 
fluctuation spectrum is therefore enhanced. Ref. [104] de- 
termined that about a third of the amplification comes 
from the increase in the radiation density while the rest 
can be attributed to the larger DM fraction. 



B. Cosmic Microwave Background 

We begin our study of the impact of the atomic DM 
scenario on the CMB by considering how the modified 
evolution of cosmological perturbations alter the spec- 
tra of temperature and polarization anisotropics. The 
CMB features caused by changes in the fluctuation evo- 
lution are more interesting than those resulting from the 
modified background cosmology since they are poten- 
tially unique to the atomic DM scenario (or, more gener- 
ally, to theories incorporating dark plasmas). To isolate 
the impact of the fluctuations, it is important to keep 
fixed quantities that only depend on the background cos- 
mology. These are z cq , 9 S , the dark energy equation of 
state parameter w, the CMB acoustic damping scale r^, 
and physical baryon density fl^h 2 . With these quantities 
fixed, we wish to determine if it is possible to distinguish 
an atomic DM model from a ACDM model containing an 
equivalent number of relativistic degrees of freedom. As 
we discuss below, the answer is positive. 

To address the issue at hand, it is instructive to first 
review the impact of relativistic neutrinos on the CMB 
in the radiation-dominated era. The key point that dis- 
tinguishes relativistic neutrinos from regular photons is 
their free- streaming nature. Indeed, while photons can 
only begin free-streaming after they decouple from the 
baryons at redshift z ~ 1100, neutrinos are generally as- 
sumed to have free-streamed since a very early epoch. 
This has two important consequences for the physics of 
the CMB [104] . First, the free-streaming of neutrinos 
causes a phase shift of the acoustic oscillations upon 
their entry into the horizon. Second, the impact of free- 
streaming radiation on the gravitational potential gener- 
ates a uniform suppression of the CMB oscillation ampli- 
tude across all multipoles. Both outcomes can be traced 
back to the facts that: (1) neutrinos propagate superson- 
ically with respect to the baryon-photon plasma and can 
thus establish metric fluctuations beyond the sound hori- 
zon of the CMB; and (2) the free propagation of neutri- 
nos sources the growth of anisotropic stress on all scales 
(including superhorizon modes) , hence affecting the grav- 
itational source terms in the photon equations of motion. 



We would like to determine how the above repercus- 
sions on the CMB change when we substitute a dark pho- 
ton for a relativistic neutrino. As with the regular pho- 
tons, the fundamental difference between the dark pho- 
tons and the neutrinos is that the former can only start 
free-streaming after they decouple from the DM (see |124j 
for a similar phenomena involving only neutrinos). This 
immediately suggests a possible way to distinguish the 
atomic DM scenario from a ACDM model containing ex- 
tra relativistic neutrinos. Indeed, if the dark photons be- 
gin free-streaming while the length scales relevant to the 
CMB arc entering the horizon, then the phase shift and 
amplitude suppression associated with the free-streaming 
of radiation will not be uniform across all the CMB mul- 
tipoles. Small-scale modes entering the horizon while the 
dark photons are still coupled to the DM will not be af- 
fected by the phase shift and amplitude suppression while 
larger scales entering after dark-photon decoupling will 
be affected, as long as they become subhorizon during 
radiation domination. These non-uniform phase shifts 
and suppressions of power across multipoles constitute 
the tell-tale signature of a relativistic degrees of freedom 
decoupling from the plasma while the Fourier modes im- 
portant to the the CMB are subhorizon. 

In Figs. [15] and |16[ we compare the temperature and 
polarization CMB power spectra for atomic DM mod- 
els having different binding energy. We also show the 
CMB spectra of a ACDM model incorporating the same 
amount of additional relativistic degrees of freedom as 
the atomic DM models such that all the spectra shown 
have the same cosmological background evolution. The 
upper panels show the spectra themselves while the lower 
panels display the relative differences between the atomic 
DM spectra and the corresponding ACDM model. We 
immediately notice that as the binding energy is in- 
creased, both temperature and polarization power spec- 
tra converge toward the ACDM model containing extra 
neutrinos. Indeed, as the binding energy of dark atoms 
is made larger, dark photons decouple earlier from the 
DM and can therefore begin free-streaming at an earlier 
epoch, making them hard to distinguish from relativistic 
neutrinos. 

For both temperature and polarization spectra, the key 
physical signatures that distinguish dark photons from 
neutrinos are clearly visible. First, let us discuss the 
amplitude suppression associated with the free-streaming 
radiation. We see that, compared with the ACDM model 
containing extra radiation, the amplitude of the CMB 
spectra for the atomic DM scenario are less suppressed, 
with the high-i multipoles being the least affected by the 
suppression. This is in line with our expectations since 
high multipoles enter the horizon before the dark pho- 
tons have the chance to significantly free-stream and are 
therefore more immune to the suppression. The atomic 
DM models with Bo — 5 keV and 10 keV clearly display 
this behavior. On the other hand, the scenarios with 
the lowest binding energies (Bp = 1 and 2 keV) exhibit 
a more complex Z-dependence when compared with the 
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FIG. 15. Comparison between the CMB angular power spec- 
trum of atomic DM and that of a ACDM model with an 
equivalent number of extra neutrinos (JV„ == 4.849 here). The 
various lines illustrate different values of the atomic binding 
energy, Bo- We fix all other dark parameters to the values in- 
dicated on the plot. The upper panel displays the TT spectra 
while the lower panel shows the fractional difference between 
the TT spectra of atomic DM and that of the ACDM model 
containing extra neutrinos. All other cosmological parameters 
are held fixed. Here, the helium fraction is fixed to Y p == 0.24. 



ACDM model. To understand this difference, we need 
to invoke the important phase shifts between the atomic 
DM and ACDM models as well as the different growth 
history of the DM fluctuations in the two scenarios. 

The effect of the phase shift is most visible for the 
Bd = 1 keV model. Indeed, since dark photons in 
this model are just beginning to free-stream when most 
Fourier modes contributing to the CMB enter the hori- 
zon, these do not experience the same phase shifts as the 
ACDM model. This can be most clearly discerned in the 
polarization spectrum (Fig. 



16). There, we see that the 



phase difference between the atomic DM models and the 
ACDM containing extra neutrinos becomes progressively 
larger toward higher multipolcs. This is exactly what we 
expect since the high multipoles enter the horizon be- 
fore they can be affected by dark-photon free-streaming, 
while smaller multipoles experience a phase shift that 
progressively converges toward the pure neutrino case as 



FIG. 16. Similar to Fig. [15] but for the CMB EE polarization 
spectra. 



I is lowered. One of the key feature of this drifting phase 
shift is that it converges toward constant values for both 
I 3> Idee and I <SC Zdeo where ^doc corresponds to the mul- 
tipole that crosses into the horizon as dark photons begin 
free-streaming. This contrasts with phase shifts caused 
by a change in the angular scale of the sound horizon 
which act multiplicatively I — > al. We show in Fig. [17] 
the damping tail of the i5i?-polarization spectrum for the 
Bp = 1 keV model. To illustrate the constant phase shift 
at I 3> ^dcc(^ 500 here), we have shifted its spectrum by 
Al ~ —8.5. We see that the phase of the shifted spectrum 
matches very well that of the ACDM model containing 
extra relativistic neutrinos, hence showing the phase shift 
does asymptote to a constant at high multipoles. 

There is an important ramification to the above ob- 
servations. Since the CMB temperature anisotropies are 
sensitive to the DM-dominated gravitational potential at 
the epoch of last scattering, any significant modification 
to the growth history of DM density perturbations will 
be reflected in the Cj T power spectrum. As we discussed 
in section |IV C[ DM fluctuations cannot grow as long as 
they are coupled to the DR. In the matter-dominated 
era, growing DM fluctuations are usually responsible for 
establishing the gravitational potentials that act as a 
restoring force to counterbalance the pressure of the pho- 
tons. If DM perturbations are prohibited to grow by their 
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FIG. 17. Illustration of the constant phase shift between an 
atomic DM model and a ACDM model containing an equiva- 
lent number of relativistic species. We show here that, in the 
atomic DM scenario, multipoles that enter the Hubble horizon 
before dark photons begin free-streaming do not experience 
the constant phase shift toward smaller / that usually char- 
acterizes models with extra relativistic degrees of freedom. 
To illustrate this, we have shifted the atomic DM spectrum 
by Al — —8.5 and observed that the phases of both spectra 
coincide. 



coupling to the DR, then the gravitational potentials can- 
not be established and the gravitational source term act- 
ing on photon fluctuations is much weaker. Observation- 
ally, this has the consequence of altering the ratios of 
the odd and even peaks of the CMB spectrum, with the 
odd peaks being suppressed on scales where DM fluctua- 
tions cannot grow. Indeed, the odd peaks correspond to 
gravity-driven compression waves which are very sensi- 
tive to the size of the DM fluctuations. The lower panel 
of Fig. [15] clearly shows the deep troughs associated with 
the damping of the odd Cj T peaks caused by the late 
kinetic decoupling of DM for the models with Bo = 1 
and 2 keV. Unfortunately, models displaying such signif- 
icant suppression of odd Cf T peaks at high multipoles 
have very low values of as and violate the Ly-cv bound 
on /cdao are therefore ruled out by observations. 

We now turn briefly our attention on how the atomic 
DM scenario affects the CMB through its impact on the 
background cosmology. Since atomic DM modifies the 
background evolution by the presence of the DR, the ef- 
fects on the CMB arc not unique but common to any 
theory incorporating extra relativistic degrees of free- 
dom. For completeness, we nevertheless discuss their 
effects here. These include a modification of the primor- 
dial helium abundance, a change to the angular scale of 
the sound horizon at decoupling, a shift of the matter- 
radiation equality epoch, an enhanced early ISW effect 
and a modified CMB damping tail. 

Primordial Helium Abundance The impact of the rela- 
tivistic dark components on the expansion rate during 
BBN tends to increase the primordial helium fraction 
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between the number of additional neutrinos at BBN 
(AN V = N v — 3) and the relativistic degrees of freedom 
of the dark sector at that epoch. Since helium rccom- 
bines before hydrogen at late times, a larger helium frac- 
tion leads to a net decrease in the free-electron fraction 
around z ~ 1100. Consequently, the photons can diffuse 
more easily out of inhomogeneities and damp tempera- 
ture and polarization anisotropies on scales smaller than 
the diffusion length. Therefore, we generically expect 
the CMB to display less power on small angular scales 
for atomic DM scenarios predicting a large primordial 
helium fraction |104l 1105) . 

Angular size of Sound Horizon, Matter-Radiation Equal- 
ity and Early ISW The presence of the dark photons af- 
fects the size of the baryon-photon sound horizon through 
its impact on the Hubble expansion rate prior to hydro- 
gen recombination. The sound horizon of the baryon- 
photon plasma is given by 
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where a* is the scale factor at recombination and c s 
is the sound speed of the baryon-photon plasma. Since 
the extra DR works to increase the Hubble rate, the net 
effect is a smaller value of r^ ^ ■ As the angular size of 
sounds horizon is given by 6 S — r^s 7 '/Z?a, where Da is 
the angular-diameter distance to the last-scattering sur- 
face, we therefore expect the CMB acoustic peaks to be 
shifted toward smaller angular scales (higher multipoles 
/). This effect is illustrated in Fig. 18 where we display 



the temperature and E-polarization spectra for different 
values of £. The shift of the acoustic peaks to smaller 
angular scales is clearly visible for both types of spectra. 
Further, we see that the temperature anisotropies are 
amplified around the first acoustic peak as £ increases. 
This is the result of the integrated Sachs- Wolfe (ISW) 
effect caused by the extra DR through its impact on the 
changing gravitational potential after recombination. In- 
deed, increasing £ brings the epoch of matter-radiation 
equality closer to that of recombination, hence increasing 
the impact of radiation on the gravitational potential at 
late times. 

Silk Damping Tail If we fix the epoch of matter-radiation 
equality, the primordial helium fraction and the angu- 
lar scale of the sound horizon at decoupling, increasing 
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FIG. 18. CMB angular power spectra in the atomic DM sce- 
nario for different values of £. We fix all other dark param- 
eters to the values indicated on the plots. The upper panel 
displays the TT spectra while the lower panel shows the EE 
polarization spectra. All other cosmological parameters are 
held fixed. Here, the helium fraction is fixed to Y p — 0.24 to 
isolate the effect from the changing sound horizon. 



the energy density of the dark photons leads to an en- 
hanced damping of the CMB anisotropics (105) 



To un- 
derstand the origin of this effect, we need to remember 
that the photon diffusion distance scales as oc H~ 5 
(see Eq. ( 63 1 ) while the angular diameter distance scales 
as H . Thus, the damping angular scale 6d = t^/Da 
effectively increases if the Hubble rate is sped up due 
to the presence of extra radiation. We therefore expect 
that as the value of £ is raised, the CMB spectrum will 
be increasingly affected by Silk damping. This effect is 
shown in Fig. [19] where we clearly observe the decline in 
amplitude associated with the increasing DR density. In 
addition, if the primordial helium fraction was allowed 
to vary according to Eq. ( 79 ), this would further increase 
the amount of Silk damping. Therefore, it is clear that 
measurements of the CMB damping tail provide strong 
constraint on £. 

In summary, beyond the impact of the atomic DM sce- 
nario on the background cosmology caused by the DR, 
we have identified four key cosmological signatures that 
distinguish the atomic DM scenario from a ACDM model 
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FIG. 19. CMB temperature power spectra in the atomic DM 
scenario for different values of £. We fix all other dark pa- 
rameters to the values indicated on the plots. We keep fixed 
throughout the redshift of matter-radiation equality and the 
angular size of the baryon-photon sound horizon at decou- 
pling. Here, the helium fraction is fixed to Y p = 0.24 to 
isolate the effect from the changing damping scale. 



containing extra relativistic neutrinos. First, the emer- 
gence of the new DAO length scale in the late-time den- 
sity field results in a minimal mass for the first DM 
protohalos that is generically larger than in the stan- 
dard WIMP paradigm. Also, as the dark photons tran- 
sition from being tightly-coupled to the dark plasma to 
a free-streaming state, they impart varying phase shifts 
and amplitude suppressions to the CMB multipoles en- 
tering the horizon. Importantly, these suppressions and 
phase shifts asymptote to constant values for I 3> ^dec 
and I <C Zdcc, a distinct feature of atomic DM that is 
not easily reproduced in the ACDM scenario. Further- 
more, we have shown that the odd Cj T peaks are sup- 
pressed on scales that enter the causal horizon before DM 
kinematically decouples. It is therefore clear that precise 
measurements on the CMB damping tail could provide 
meaningful constraints on the parameter space of atomic 
DM. We should however keep in mind that the modi- 
fied evolution of DM and DR fluctuations can only affect 
the CMB if the dark sector kinetic decoupling happens 
close enough to the epoch of last scattering. As such, a 
non-detection of these signatures effectively puts a lower 
bound on the redshift of kinetic decoupling which itself 
depends on a combination of ao, Bp, mp, and £. 



VI. ASTROPHYSICAL CONSTRAINTS ON 
ATOMIC DARK MATTER 

As the Universe expands and cools down, non-linear 
structures begin to emerge and eventually form present- 
day astrophysical objects such as galaxies and clusters of 
galaxies. The internal dynamics of these objects is deeply 
influenced by the microphysics governing DM because the 
latter contributes the vast majority of the mass inside 
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these objects. Since the atomic DM scenario naturally 
incorporates new interactions in the dark sector, it is 
important to discuss its implications for the dynamics of 
DM halos on a wide range of scales. 

On one side, observations indicate that the DM halos 
of elliptical galaxies and clusters display a triaxial ellip- 
soidal shape [34j I125H133] . This indicates that the relax- 
ation (thermalization) time of their DM halos is much 
longer than the typical dynamical time of these celestial 
objects. From this observation, one can obtain a bound 
on the elastic collisional rate of DM particles. Similarly, 
detailed studies of the Bullet Cluster provide direct con- 
straint on the self-interaction cross section of DM |134j . 
However, these constraints have been shown to be weaker 
than those derived from the halo ellipticity |63j and we 
therefore do not consider them here. 

On the other side, due to the rich internal structure of 
dark atoms, the atomic DM scenario inherently includes 
dissipation mechanisms that could potentially have dra- 
matic effects on the dynamics of DM halos. Indeed, 
collisional excitations of dark atoms followed by dark- 
photon emissions provide a cooling mechanisms for DM 
that could drastically alter the internal structure of halos. 
Further, the chemistry of the dark sector naturally pro- 
vides other heat-dissipation mechanisms such as molecu- 
lar cooling. Fortunately, demanding that DM particles be 
effectively collision-less to preserve the structure of halos 
strongly hinder the efficiency of the atomic dissipative 
processes, rendering them mostly irrelevant. Similarly, 
the requirement that very few actual particle collisions 
take place in the dark sector most likely shut off any 
chemical reactions inside halos. We therefore do not fur- 
ther consider the possibility of dark chemistry. 



A. Ellipticity of DM Halos 

Collisions among the different constituents of the dark 
sector lead to a redistribution of the linear and angu- 
lar momentum inside a DM halo. This tends to erase 
velocity correlations on halo scales and brings the DM 
halo closer to isothermality, effectively altering the halo 
structure and making it roughly isotropic. Since observa- 
tions indicate that halo iso-density contours are elliptical 
and thus significantly deviate from isotropy, we conclude 
that the DM thermalization time is much longer than the 
dynamical timescale of halos. This in turns tightly con- 
strains the rate of particle collisions in which a significant 
amount of momentum is exchanged. 

In the atomic DM scenario, we need to consider col- 
lisions between the dark atoms, the dark electrons, and 
the dark protons. In total, this amounts to six differ- 
ent types of collisions: H-H, H-e, H-p, e-e, e-p, and 
p-p. Here, H stands for the dark atoms. Since we are 
most interested in dark sector that are mostly neutral 
when non-linear structures form, the contribution from 
the first three types of collision is expected to dominate. 
To compute the rate of collisions with large momentum 



exchange, we need the momentum-transfer cross section 
which is defined by 
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where da/dtt is the elastic differential cross section. 
A detailed computation of the H-H, H-e and H-p 
momentum-transfer cross section is beyond the scope of 
this paper, but we can use the extensive literature (see 
e.g. |135j ) on the corresponding visible-sector cross sec- 
tions to derive some general properties. We refer the 
reader to appendix [D] for further details. To a good 
approximation, the H-H, H-p, and H-e momentum- 
transfer cross sections are 
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where [iDe = rn£,m e / {mo + m e ) ~ m e and v is the rel- 
ative velocity of the colliding dark ions or atoms. The 
ion-ion differential cross sections are simply given by the 
Rutherford scattering cross section. Performing the inte- 



gral in Eq. 82 leads to 
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where 0~^ n is the minimum scattering angle. We take 
this angle to be given by the Debye screening length Aoe 
of the dark plasma through the relation 
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and where n e is the number density of free dark elec- 
trons inside a DM halo. The momentum-loss rate of dark 
species i upon collisions with species j in a DM halo is 
given by 
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where n- } = nj (r) is the number density of specie j inside 
a halo, v is the relative velocity of the collision, Api is the 
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momentum loss in a single collision and f(v) is the veloc- 
ity distribution which we take to be locally Maxwellian 
with velocity dispersion vq 
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The momentum lost by particle i upon colliding with 
particle j is 
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where pi is the momentum of particle i. Normally, the 
momentum loss is weighted by a factor (1 — cos 6), but 
we have absorbed this factor into the momentum-transfer 
cross section. The rate of momentum-changing collisions 
between dark species i and j is defined as 



= Pr 
Pi 
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where pi is the momentum of particle i averaged over the 
velocity distribution given in Eq. ( 92 ) . With the help of 



the above momentum-transfer cross sections, these colli- 
sions rates are 
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In the above, T(x) is the Gamma function. The rates 
for the opposite processes are obtained by rescaling the 
above rates with the appropriate mass densities 
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In the case for which m e <C to p , ttid, we see that 
Tpe *C r ep and Tn e <C Tne if the abundance of dark 
atoms and dark ions is roughly similar. This is reason- 
able: a heavy dark proton needs to scatter off many light 
dark electrons before its momentum is significantly af- 
fected. Conversely, a dark electron's momentum can be 
dramatically changed by a single collision with a dark 
proton or a dark atom. However, since the dark electrons 
are usually much lighter than the dark protons and there- 
fore carry a small fraction of the overall halo mass, these 
high-momcntum-transfcr collisions do not affect the el- 
lipticity of DM halos. Therefore, we do not consider any 
further the scattering of dark electrons (that is, we ne- 
glect ee, ep and eH collisions). In the case where the 
dark proton and the dark electron are nearly degenerate 
in mass, constraints on dark proton collisions naturally 
engulf the dark electron constraints so we can still neglect 
them. 

Detailed simulations of self-interacting DM halos have 
shown that DM particles forming the bulk of the matter 
density can undergo up to 10 hard scatters in a Hubble 
time before the ellipticity of halos is adversely affected 
[IM1CE3I]- Therefore, the ellipticity of DM halos is pre- 
served if we have 
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where r co n is the overall hard-scatter rate of the dark 
sector. The total collision rate inside a halo is the sum of 
the individual rates weighted by the relative abundance 
of each dark species 
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To compute the above constraint, we need to specify the 
density profile of both dark ions and dark atoms as well 
as the velocity dispersion vq (r) . The main difficulty here 
is that the actual density profile of each dark constituents 
is itself determined by their collisional rates. For exam- 
ple, if the dark ions undergo many hard scatters among 
themselves within the typical dynamical time of a galaxy, 
they settle into their own isothermal density profile [43], 
while the neutral dark atoms maintain their CDM-like 
profile. This implies that the local ionized fraction inside 
a DM halo can vastly differs from the background values 
computed in section IIIB| (by several order of magnitudes 
in some cases). To complicate matters further, collisions 
between dark atoms in the central region of halos likely 
lead to the formation of a cored density profile. There- 
fore, it is clear that detailed iV-body simulations would 
be required to determine the actual density profile of each 
dark constituent. 
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We can nevertheless derive conservative constraints on 
atomic DM by making some simple assumptions. We 
first assume halos to be locally neutral such that n e (r) = 
n p (r). We further assume the ionized component of the 
DM follows the total density profile, that is, n e (r) = 
n p (r) = x D n D (r). Similarly, n H (r) = (1 - x D )n D (r). 
We take the halo ionized fraction (xd) to be equal to 
the late-time ionized fraction of the background, that is 
xd = %d{z = 0). It is important to emphasize that 
the above prescription is not accurate since the ionized 
and neutral components generally have different density 
profiles [43] and the ionized fraction inside halos may sig- 
nificantly differ from the overall background cosmological 
values. As such, the bounds derived from these assump- 
tions should therefore be considered as upper limits on 
how limiting the ellipticity constraint can be. The actual 
bounds on the parameter space of atomic DM are poten- 
tially much less constraining. Nevertheless, the above 
treatment is expected to be fairly accurate for ij; « 1 
and Xd ~ 1. 

While it is true that the most stringent constraints on 
halo ellipticity (see e.g. [63 ) come from the inner part of 
galactic halos, it is unclear how these bounds apply to 
atomic DM. Indeed, ellipticity constraints strongly de- 
pend on the local DM density in the inner central re- 
gion of halos. The latter is usually obtained by fitting 
a predefined halo profile (e.g. NFW) to the data |138j . 
Since atomic DM generally admits a different halo shape, 
we cannot blindly apply these results to the atomic DM 
scenario. We therefore resort to typical values of the 
DM density and velocity dispersion inside galactic ha- 
los. Explicitly, we evaluate the above constraints using 
pn ~ 3, 1, and 0.3 GeV/cm -3 . We take the velocity dis- 
persion to be Vq ~ 250 km/s. The number density of DM 
particles is then given by jid — Pd/ht-d- 

We display in Fig. [20]the constraints oaaj) and mr> 
for six different values of the atomic binding energy. 
We display the three disfavored contours corresponding 
po — 3,1, and 0.3 GeV/cm~ 3 . We observe that for 
lOkeV < Br> < 100 keV, there is little parameter space 
for which the dark sector is mostly neutral and is not 
in tension with the ellipticity constraint. As the bind- 
ing energy is increased above this threshold, a large al- 
lowed parameter space opens up since the atomic geomet- 
ric cross section rapidly decreases as the binding energy 
climbs up in value. Below 10 keV, another unconstrained 
region opens up at large coupling constant values and 
masses. This is caused by the collision energy approach- 
ing the excitation threshold of the dark atoms. As this 
limit is approached, the momentum-transfer cross sec- 
tion for atom-atom scattering become more and more 
suppressed, while the inelastic channels start to growth 
in importance. The inelastic cross sections are however 
much smaller than the elastic ones (~ ira^ D instead of 
~ 10 2 7rciQ D ) and we thus do not expect these collisions 
to severely affect the ellipticity of haloes. 

Overall, the ellipticity constraint unambiguously disfa- 
vors dark atoms with mass < 1 TeV, unless Bp > 1 



MeV. We however reiterate that the constraints shown in 



Fig. 20 are conservative, especially in the regions where 
1% ^5 xd ^5 99% for which detailed simulations of dark- 
atomic halos will likely be required to assess the valid- 
ity of the bounds. Indeed, if the typical hard-scatter 
timescale of dark ions is shorter than the dynamical time 
of a galaxy we can safely assume that the dark ions set- 
tle into a separate isothermal density profile [33] . In the 
presence of a sizable population of neutral dark atoms, 
the overall DM halo can however still display significant 
ellipticity if the dark atoms are themselves mostly colli- 
sionless. Since the isothermal ionic halo is typically much 
more diffuse |43j than the neutral one, ion-ion and ion- 
atom interactions are probably much more suppressed 
than our naive estimate suggests. This is likely to open 
up much of the parameter space for 1% < xd < 99%. 

In the case where the dark ions form the majority of 
the matter density inside halos, the ellipticity constraints 
derived in the context of hidden charged DM apply [62 \ 
IS3"] . In this scenario, dark ions cannot settle into an 
isothermal profile without utterly violating the ellipticity 
of halos. Viable models can nevertheless be constructed 
by suppressing the dark fine-structure constant and by 
considering very large ions masses (such that the number 
density is low). Fortunately, as ao is increased and m_o 
is decreased, more and more ions are bound into dark 
atoms hence softening the ion-scattering bound. In fact, 
since the early recombination rate scales as ~ a D while 
the ion- ion hard-scatter rates scale as ~ a D , the ionized 
fraction drops faster than the ion-ion momentum-transfer 
rate is increasing as the the dark fine-structure constant 
is dialed up. Thus, the atomic DM scenario naturally 
provides a way to evade the halo ellipticity constraints 
on hidden-charged DM. 

In summary, while detailed simulations of DM halo 
formation in the atomic DM scenario are likely required 
to determine the exact constraints, it is clear that dark 
atoms lighter than 1 TeV are likely to lead to collisional 
DM that is in tension with the observed internal structure 
of halos, unless Bjj > 1 MeV. 



B. Cooling of DM Halos 

Collisional excitation of dark atoms followed by the 
emission of a dark photon provides a natural cooling 
mechanism for DM. Since DM is observed to have very 
different properties than baryons (which are allowed to 
cool), collisional cooling of dark atoms must be sup- 
pressed. The easiest channel to excite a dark atom is 
through the hyperfine transition. The hyperfine splitting 
is given by 
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where g e ~ 2 and g p ~ 2 are the gyromagnetic ratio 
of the dark electron and dark proton, respectively. This 
transition can be excited when a dark atom in a spin 
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FIG. 20. Halo ellipticity and Ly-a forest constraints on the parameter space of atomic DM. We display the constraint for two 
values of the atomic binding energy. The allowed region (orange) is clearly indicated on the plots. The blue, yellow and purple 
contours (outermost to innermost) display the disfavored regions when po — 3, 1, and 0.3 GeV/cm -3 , respectively. The light 
and dark gray regions correspond to minimal halo masses M m i n > 9.3 x 10 10 Mq and M m i n > 10 x 10 13 Mq, respectively. 
The white dashed lines indicate contours of constant minimal halo mass. We also show the contours (solid black) of constant 
background ionization fraction for xd = 1% and xo = 99%. 
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singlet state collides with a dark ion or atom and under- 
goes a spin-flip to the triplet state. The cross section for 
this process is somewhat smaller than the elastic scatter- 
ing cross section [135] . Therefore, it is most likely that 
all regions of parameters space where hyperfine emission 
leads to significant cooling are already ruled by the ellip- 
ticity constraint. To verify this, we compute the typical 
timescale for a dark atom to lose an 0(1) fraction of its 
kinetic energy due to hyperfine emission. Assuming that 
this process is approximately governed by the geometric 
cross section, this timescale is 
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Demanding that this timescale be longer than the age of 
the Universe leads to no new constraints beyond those 
already plotted in Fig. |20] 

Beyond the hyperfine transition, the other dissipative 
process that is relevant for dark atoms is the Is — > 2s, 2p 
collisional excitation. Heuristically, such inelastic colli- 
sions can only happen when the timescale of the collision 
is shorter than the typical timescale of the dark atoms 
(i.e. adiabaticity is violated). The typical velocity of a 
dark electron inside an atom is v e ~ an, leading to an 
atomic timescale td ~ &o,£>A?e- On the other hand, the 
collision timescale is of the order t co u ~ o-o,d/ v , where v 
is the relative velocity of the collision. Taking r co n < fjj , 
inelastic collisions are only possible when 
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According to Fig. 20 all allowed regions of parameter 
space displaying a large fraction of neutral dark atoms 
have Od > 10~ 2 . Taking v to be approximately equal to 
the typical velocity dispersion inside a DM halo, we have 
v ~ 10~ 3 for galactic halos while v ~ 10~ 2 for galaxy 
clusters. It is thus clear that inelastic collisions can play 
no major role inside galactic halos while they could play 
a marginal role inside clusters. In the latter case how- 
ever, the small number density of dark atoms coupled 
with the typically small inelastic cross sections (see e.g. 
[139] ) likely render inelastic collisions completely negligi- 
ble inside clusters. We therefore conclude that inelastic 
processes do not constrain the atomic DM scenario be- 
yond the regions of parameter space already ruled out by 
the cllipticity and the matter power spectrum bounds. 



VII. DIRECT DETECTION 

Most previous works on atomic DM [4"2Tl4"4l 1140] have 
focused their attention on the potential direct-detection 
signatures. Since atomic DM can naturally scatter in- 
elastically, it offers a mechanism to reconcile the annual 



modulation seen by DAMA [141] and CoGeNT [142] with 
the null signal of the CDMS [1431 QH] and XenonlO 
[145H147] experiments. The actual direct-detection re- 
sults strongly depend on how atomic DM couples to the 



SM. In this section, we compare the atomic DM mod- 
els which provide a good fit to the direct-detection data 
to our cosmological and astrophysical constraints derived 
above. 

In Ref. [42], the authors considered a model in which 
dark fermions are axially coupled to a broken U(l) which 
mixes with the SM hypercharge. They found that dark 
atoms having an hyperfine splitting E^f ~ 100 keV, 



70 < 



< 200 GeV and 2.2 > 



> 1.6 GeV can 



provide a good fit to the modulated spectrum of DAMA 
while evading the constraints from other direct-detection 
experiments. This corresponds to an atomic DM sce- 
nario with 27 < B D < 45 MeV, 0.16 < a D < 0.24 
and 72 < vtld < 201 GeV. Extrapolating the constraints 



shown in Fig. 20 to the appropriate binding energies, we 
find that these values lie on the edge of the ellipticity 
bound and are therefore marginally allowed by our con- 
straints. 

In Ref. |43j . the authors considered a similar broken 
axial U(l) model which mixes with the SM hypercharge. 
In this work, they show that an atomic DM scenario with 
Ehf ~ 5 — 15 keV with dark fermion masses m e ~ m p ~ 5 
GeV could reconcile the CoGeNT data with the count 
rate seen by the CRESST experiment [148] . This cor- 
responds to an atomic DM scenario with an ~ 0.03, 
Bd ~ 1 MeV and mn ~ 10 GeV. This model is in seri- 
ous tension with the halo constraint shown in Fig. [20] and 
is therefore likely to lead to collisional DM which would 
dramatically affect the structure of halos. 



In Refs. [HI I140j . the authors considered atomic DM 
models where the dark photon mixes with the regular 
photon. This effectively gives small standard electro- 
magnetic charges to the dark fermions. The authors con- 
cluded that scenarios with an ~ 0.04 — 0.06, Bd ~ 2 — 5 
MeV and ttid ^6—10 GeV predict the right cross sec- 
tion to explain the DM events recorded by the CoGeNT 
collaboration. These values are however clearly in ten- 
sion with the halo ellipticity bound. Nevertheless, since 
they lie close to the boundary with the allowed region 
and since our constraints are conservative, detailed N- 
body simulations will be required to assess whether this 
model is ruled out or not. 

We observe that dark-atom models capable of explain- 
ing the current direct-detection experiments are at best 
marginally allowed by our conservative astrophysical con- 
straints. This stems from the fact that direct-detection 
experiments generally favor vtld ~ 1 — 10 GeV while the 
halo constraint typically imposes > 1 TeV, unless 
Bd > 1 MeV. From the point of view of direct detection, 
the main challenge for building a successful atomic DM 
theory is to avoid strong elastic scattering off nuclei while 
allowing for enough inelastic collisions to explain both 
DAMA and CoGeNT. The study of atomic DM models 
that could be in agreement with both the direct-detection 
data and the astrophysical constraints are left for future 
work. 
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VIII. DISCUSSION 

We have presented in this work a thorough study of 
a dark sector made up of atom-like bound states. This 
model naturally incorporates extra relativistic degrees of 
freedom whose presence are currently favored by CMB 
experiments. The observed primordial abundances of 
light elements constrain the temperature of DR to be 
somewhat cooler than the CMB, with an upper bound 
given approximately by Td/Tcmb < 0.8 at z = 0. We 
have revisited the atomic physics necessary to describe 
the processes of dark recombination, thermal decoupling, 
and kinetic decoupling. We find that in some cases, the 
inclusion of physics beyond that of standard atomic hy- 
drogen is required to properly describe these transitions. 
In particular, we find that bound-free processes such as 
photo-ionization heating and photo-recombination cool- 
ing are key to determine the thermal decoupling epoch 
of weakly-coupled atomic models. For strongly-coupled 
dark atoms, we have shown that the addition of Rayleigh 
scattering can significantly delay the kinetic and thermal 
decoupling of DM. 

We have solved the linear cosmological perturbation 
equations, taking into account the interaction between 
DM and DR and showed that atomic DM can go through 
various regimes as time evolves and as the dark param- 
eters are varied. In particular, the DR pressure leads to 
strong acoustic oscillations for Fourier modes entering the 
Hubble horizon prior to the dark drag epoch. Further, 
diffusion and acoustic damping severely suppress the am- 
plitude of DM fluctuations on scales shorter than the 
sound horizon at kinetic decoupling. At late times, these 
features remain imprinted on the matter power spectrum. 
Importantly, the atomic DM scenario introduces the new 
DAO length scale in the density field which basically de- 
termines the minimal DM halo mass. We have shown 
that observations of Ly-a forest flux power spectrum put 
an upper bound on the size of the DAO scale and rule 
out a large fraction of atomic DM models with Bo < 1 
keV. 

We have performed a detailed study of the impact of 
the atomic DM scenario on the CMB. We have deter- 
mined that the largest impact on the CMB in this model 
is due to the presence of dark photons. If dark photons 
begin free-streaming at a very early epoch, their impact 
on the CMB is likely indistinguishable from that of extra 
relativistic neutrinos. On the other hand, if dark photons 
decouple from DM when Fourier modes relevant for the 
CMB are inside the horizon, then the atomic DM scenario 
predicts CMB signatures that are difficult to reproduce 
with only relativistic neutrinos. These signatures include 
non-uniform phase shifts and amplitude suppressions of 
the temperature anisotropy power spectrum. The clear- 
est CMB signature of the atomic DM scenario is a non- 
uniform phase shift of the polarization power spectrum 
that asymptotes to a constant at both large and small 
scales. 

These signatures can only be present if the dark pho- 



tons form a sizable fraction of the radiation energy den- 
sity (£ ^ 0-6) and if they decouple at a late enough red- 
shift, which usually requires Bjj < 10 kcV. These rela- 
tively low atomic binding energies are however strongly 
constrained by both the Ly-a data and the halo ellip- 
ticity requirement. For the remaining allowed parame- 
ter space, dark photons generically decouple too early to 
have an impact on the CMB that is noticeably different 
from that of standard relativistic neutrinos. It is there- 
fore unlikely that the CMB can be used to learn about the 
simplest atomic DM model discussed here. Nevertheless, 
in a more general model where dark radiation only cou- 
ples to a certain fraction of the DM, it is possible that 
the former could significantly affect the CMB without 
modifying the small-scale matter power spectrum. This 
is an intriguing possibility that we leave for future work. 

The strongest constraint on atomic DM comes from re- 
quiring that DM is effectively collision-less inside galac- 
tic halos. This stems from the fact that atoms naturally 
have much larger geometric cross sections than point par- 
ticles. Since this geometric cross section scales as Bp 2 , 
the halo constraint favors a large binding energy, with 
models having Bu > 10 MeV largely unconstrained. Un- 
surprisingly, this corresponds to the regime where atomic 
DM closely resembles a standard WIMP particle. For 
Bjj < 1 keV, the ellipticity and Ly-a constraints favor a 
dark sector that is mostly ionized at all times. Our con- 
straints in these regions of parameter space are similar to 
those found in Refs. [621163] - Importantly, we have shown 
that inelastic processes are inefficient at dissipating the 
energy of dark atoms in all the regions that are not al- 
ready ruled out by the halo and the small-scale power 
spectrum bounds. Overall, it is clear that preserving the 
observed ellipticity of halos severely limits the possible 
interactions between DM particles. 

Since the atom-atom, ion-atom and ion-ion cross sec- 
tions are all velocity dependent, atomic DM can provide 
enough interactions to smooth out the central region of 
small satellites galaxies while retaining the ellipticity of 
large galactic and cluster halos. Furthermore, atomic DM 
scenarios predicting a minimal halo mass in the range 
M min ~ 10 8 - 10 10 M have the potential to affect the 
faint-end of the galaxy luminosity function and bring it 
in line with observations. In addition, the presence of an 
ionized component naturally makes the halo more diffuse 
and could potentially alleviate the so-called "too big to 
fail" problem (THE)]. Interestingly, an atomic DM model 
with B D ~ 5 keV, m D ~ 80 TeV and a D ~ 0.02 could 
possibly address all three problems affecting dwarf galax- 
ies. Indeed, this model has a late-time global ionized frac- 
tion of xd — 0.2, a minimal halo mass of M m - m ~ 8 x 10 8 
M Q and lies at the boundary of the collisional constraint, 
meaning it could contain enough interaction to form cores 
in galactic halos while retaining their overall ellipticity. 
A detailed numerical study will be necessary to assess 
the success of this model. 

One issue that we have not touched upon in this work 
is the possibility of dark magnetic fields. These have the 
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potential to significantly alter structure formation unless 
there is a mechanism that naturally suppresses their am- 
plitude and range. The simplest way to achieve such 
suppression is to break the U(1)d gauge force by intro- 
ducing a small dark-photon mass. This is an intriguing 
possibility since, although a minute mass is needed to 
quench large-scale magnetic fields, a mass whose Yukawa 
length is comparable to the size of a dark atom could 
considerably alleviate the halo constraint by shrinking 
the size of the atomic cross section and by limiting the 
range of ion-ion interactions. One would however have 
to revisit the atomic physics described in section |III| to 
make quantitative predictions about the thermal history 
of such a dark sector. We leave the study of such a model 
to future work. 

Direct-detection data favor light atomic DM candi- 
dates which are in tension with the halo constraint. Nev- 
ertheless, the direct-detection signatures strongly depend 
on how dark atoms couple to the SM. Whether a success- 
ful model which agrees with both the halo and direct- 
detection constraint could be constructed remains an 
open question. While an explicit model would have to 
be specified in order to make quantitative predictions, 
we emphasize that dark atoms are still a viable DM can- 
didate from the point of view of both astrophysics and 
direct detection. 



Appendix B: Effective Number of Relativistic 
Degrees of Freedom 

The effective number of degrees of freedom contribut- 
ing to the entropy density of the DS at a temperature Tp 
is 
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where y e = ^ and the first term corresponds to the 
dark-photon contribution while the term in square brack- 
ets denote the dark-electron contribution. Performing 
the integrals, we obtain, to a very good approximation, 



g* s>D (T D ) ^2 + ^(1 



Ve 



1.394\0. 247 -0.277^ 



(B2) 



Appendix C: Thermal Rates 



In this appendix, we compute the rates governing the 
energy exchange between the DM and the DR. 
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1. Photo- Ionization Heating Rate 



Appendix A: Relations between Dark Parameters 

In this appendix, we relate the dark electron mass m e , 
the dark proton mass m p , and the atomic reduced mass 
fj,D to the dark sector parameters. 
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The photo-ionization heating rate is given by 

<' Lnl ° D L 



n p _i = V n„iB D i 
„ / Jo 



e B D (K 2 +n- 2 )/T D _ 1 



d( K 2 ), (CI) 



where the sum runs over all atomic states specified by the 
quantum numbers n > 2 and I < n— 1, n n i represents the 
number density of dark atoms in state nl and 7 n ;(/t) is 
defined in Eq. (17 1. We do not include bound- free tran- 



sitions to and from the ground state since recombination 
directly to the Is state results in the emission of an ener- 
getic photon that is immediately reabsorbed by a nearby 
neutral atom, hence resulting in no net cooling or heat- 
ing. We assume that the occupations numbers of excited 
states are in Boltzmann equilibrium with the 2s state 



(C2) 



where g n i is the degeneracy of the energy level with quan- 
tum numbers nl. This approximation is valid as long as 
Bd{Td) ^ H. Performing the momentum integration in 
Eq. (CI) and computing the sum over atomic states up 
to 7i max = 250 yield a photo-ionization heating rate of 
the form 

q,3 J>2 b d 

n p _i(T D ) = -^x 2s n D e *n>F p -i(T D /B D ), (C3) 

-F p _i(y) is a dimensionless universal function encoding 
(A2) the remaining of the temperature dependence of the 
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photo-ionization heating rate. It is shown in Fig. [21] 
This function is well-fitted by 



F P -i(y) 

for 4 x 1CT 4 < y < 10 2 . 
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(2.012 + y 0.2412)6.55' 



(C4) 



10.0 

'q 5.0 
cq 

^ 1.0 

_ o. 

ki 0.5 



0.1 



10" 4 



0.001 



0.01 0.1 

T D /B D 



10 



100 



FIG. 21. Universal dimensionless fitting functions F p _i for 
the photo-ionization heating rate, plotted as a function of 
To/Bd- 
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FIG. 22. Universal dimensionless fitting function F p _ r for 
the photo-recombination cooling rate, plotted as a function 
of T D /B D for 4 different constant values of T DM /T D . 



thermal dark electron obeying a Maxwell-Boltzmann ve- 
locity distribution. The rate for thermal Bremsstrahlung 
cooling is [HOI CEEE] 



n 



ff-c 



16a 3 D g s ^/2TrT DM x 2 D n 2 D 
(3 M Z3) 3 / 2 
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2. Photo-Recombination Cooling Rate 



The rate of photo-recombination cooling is given by 



Tlp-r = 



(2tt) 3 / 2 B d 



-D' l D 



n , / 



(C5) 



x / e 



-Bdk /Tdm k 2 



k 7„/(k) 



x [l + ZBB^^+n- 2 ),^)] d( K 2 ). 

As before, we do not consider recombination directly to 
the ground state. Computing the momentum integral 
and summing over all atomic states up to n max = 250 
yield a photo-recombination cooling rate of the form 



n p _ r 



2a 3 D V2vrT J 
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Bd ' Td 
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where F p _ r is a dimensionless universal function. We 



illustrate its behavior in Fig. [22} 



3. Free-Free Cooling and Heating 

The free-free emission (absorption) process refers to 
the emission (absorption) of a dark photon by a dark elec- 
tron due to its acceleration in the Coulomb field created 
by a dark proton. For atomic DM, we are mostly inter- 
ested in the emission and absorption of dark radiation by 



where gs is the free-free velocity-averaged and frequency- 
averaged Gaunt factor, which is of order unity. The in- 
verse process of thermal free- free absorption is character- 
ized by the following absorption coefficient [151] 

o 3 2 2 CI _ -2-kv/T dm \ K 
„tt(rp \ za D n D x D\ L e )9v /na \ 

(^>\xdY ii \/^Tdmv a 

where g^ is the velocity-averaged free-free Gaunt factor 
and v is the frequency of the incoming dark photons. 
The total amount of energy absorbed through free-free 
interactions per unit volume per unit time is then 
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where B v {Td) is the Planck function. To evaluate this 
integral, we first note that the integrand only depends on 
Tdm itself and on the fractional temperature difference 
between the DM and DR, e = (T D — T DM )/T D 



,2-Kv/Tnu _ l 



,-2-ku/Tdm _ p -i™/T D M 



e 2-Kv/T D _ i e -2wue/T D M — e -2vv/T DU ' 

(CIO) 

In the limit of quasi thermal equilibrium, e -C 1, we can 
expand the integrand as a power series in e and com- 
pute the frequency integral order-by-order in the frac- 
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tional temperature difference. This leads to 



Appendix D: Atomic Cross Sections 
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(Cll) 



12tt 



(6 + ir 2 e{l + 2e) - 6e 2 C(3) + 0(e d )), 



where C( x ) is the Riemann zeta function. The net rate 
at which DM gains energy due to free-free interactions, 
n ff ee n ff _ h - riff-c, is the given by 



n ff 



16a 3 D gtf^2TTT DM x 2 D n 2 D 
(3/x B ) 3 / 2 
Tr 2 e(l + 2e) -6C(3)e 2 
6 



For numerical computation, we take gs ~ 1.3. 



(C12) 



4. Rayleigh Heating 

After the onset of recombination, dark photons can 
transfer energy to the DM sector via Rayleigh scattering 
off neutral dark atoms. The scattering cross-section for 
this process is given by 



(C13) 



a R {u) rj a T:D — 



where a-r,D = 87ra|,/(3mg) is the dark Thomson cross- 
section and fhya is the frequency of dark Lyman-a pho- 
tons. This expression is valid for v -c Vhya, which is real- 
ized for the large majority of dark photons after the onset 
of dark recombination. The net rate of energy transfer 
between dark photons and DM due to Rayleigh scatter- 
ing is 

a K {u)B y {T D )AE K ^-, (C14) 

where -B„(Td) is the Planck function and AE R is the 
average net energy gained by a dark atoms in a Rayleigh 
scattering event. It is given by 
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(2m> - wT DM ) 



(C15) 



where w = 75600^(9) /ir s is a constant that can be de- 
termined using detailed balance. Substituting Eq. (C15) 



into Eq. ( C14 ) and performing the integral over frequency 



yields the effective heating rate 
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where £(y) stands for the Riemann Zeta function. This 
expression is valid for Tjj <C Bjj which is always realized 
when a large population of neutral dark atoms is present. 



In this appendix, we justify the form of the cross sec- 
tions used to describe collisions of DM particles. 



1. H-H Cross Section 

We consider the elastic scattering of two identical 
hydrogen-like atoms. This problem requires one to solve 
for the joint wavefunction of the two atoms interacting 
via the singlet gerade and triplet ungerade molecular po- 
tentials |135l 1153] . The elastic differential cross section 
for indistinguishable hydrogen-like atom is given by 
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where the subscripts "s" and "t" refer to the singlet and 
triplet states, respectively. The singlet and triplet con- 
tributions take the form 
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where the uppermost sign is for the singlet and the lower 
sign for the triplet. Here, f s ,t(6) is the scattering ampli- 
tude for the centre-of-mass (CM) scattering angle 9 and 
is given by 



/ s ,tW = 2^E( 2Z+1 )( e 



2iSf J 



l)P;(cos( 



(D3) 



where k is the wavenumber in the CM frame, / is the 
orbital angular momentum, P;(cos6*) is the I th Legendre 
polynomial, and is the scattering phase shift from 
either the singlet or the triplet molecular potential. From 
Eqs. (D2) and (D3), one can obtain the important result 
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(D4) 



are the momentum-transfer and 
the total elastic cross section, respectively. This result 
can easily be proven using the orthogonality of the Leg- 
endre polynomial and follows from the symmetry of scat- 
tering two identical atoms. 

All the information about the atomic physics is en- 
coded in the phase shifts 5^. These can be found by 
solving the Schrodinger equation 



1(1 + 1) 



+ 2 f i[E-v a ' t (R)] 



Ru*'\R) = 0, 



(D5) 



dR 2 R 2 
such that the asymptotic behavior at large R is 

u{\R) ~ R- 1 sin (kR - Ztt/2 + 5^). (D6) 
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Here, \x is the reduced mass of the colliding atoms, 
V s,t (R) is the interatomic potential in the singlet or 
triplet state, E is the relative energy of the collision and 
k = \j2[iE. In the semi-classical limit, the phase shift 
can be written as [1351 1154] 



JfS.t 



fc 2 - 2^y s - t (i?) - 



R 2 



dR 



1(1 + 1) 
R 2 



dR, 



(D7) 



where i?o and R' are the outermost zeros of the respec- 
tive integrands. 

Now focusing on the atomic DM scenario, we would 
like to understand how the phase shifts are affected when 
the dark parameters ajj, Bp and mp are varied. These 
parameters enter the phase shifts through their contribu- 
tions to the wave number k and to the molecular poten- 
tials V s,t (R). While it is out of the scope of this appendix 
to compute the molecular potentials of dark atoms, we 
can nonetheless extract general properties of these po- 
tentials by studying the case of standard molecular hy- 
drogen. To gain insights into the scaling of the molecular 
potential with the dark parameters, let us briefly consider 
the Lennard- Jones potential [155j 



only depend on two dimensionless combinations of the 
dark parameters 
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where we used Eq. (|D9|) in the last line. Since Vh/-Bh 



is known, once could then use Eq. (D10) to obtain the 



phase shifts for dark atoms. However, Eq. (Dll ) suggests 



that one could obtain approximate <5; ; £> directly from the 
baryonic hydrogen phase shifts by rescaling the collision 
energy as 
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up 
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Of course, this is only an approximation since the molec- 
ular poten tial would also have to be rescaled according to 
Eq. (Dill. Moreover, the molecular potential itself de- 



pends non-trivially on ap and \ip and would have to be 
modified accordingly. Nevertheless, studies of deuterium, 
tritium |156j and positronium |157j self-scattering sup- 
port the above approximate rescaling. We therefore take 
our atom-atom momentum-transfer cross section to be 



V(R) = e 
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~R 
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where the first term accounts for the short-range repul- 
sive force due to Pauli blocking while the second term de- 
scribes the long-range behavior of the interaction. Here, 
e is a constant setting the depth of the potential well and 
R m is the distance where the potential admits a mini- 
mum. 

On purely dimensional ground, the typical energy scale 
of the molecular potential should be proportional to the 
atomic binding energy, that is, e cx Bp. Also, the typical 
atomic separation R m should be of the order of the dark 
atomic Bohr radius, R m oc ao,_D = Zap/Bp. We can 
thus write 



V D (R/a 0tD ) 
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(D9) 



where the "H" subscript stands for the standard bary- 
onic hydrogen. Defining the dimensionless distance r = 
R/ a o,D, we can write Eq. (D7) as 



J LP 
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dr 
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where we have substituted the two-atom reduced mass 
fi = mo/2. We therefore observe that the phase shifts 



H H 



E] 



a 0,P 

a 



(D13) 



x a 



H-H 



ap 



mp 
m H 



Bh 
Bp 



E 



where the leading factor accounts for the change in the 
geometric cross section. A numerical fit to the calcula- 
tions found in Refs. [1351 1156j leads to 



a^ H (E CM )-lO z 7Tai 



(D14) 



which is valid for 10 -3 < E C m/B h < 10. At high en- 
ergy, inelastic collisions begin to dominate over the elastic 
channel and this explain the much steeper energy depen- 
dence in that regime. Using Eq. (D13l, we can finally 



obtain an approximate momentum-transfer cross section 
for the neutral dark atoms 
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where we have use the definition of binding energy to 
simplify the energy rescaling. Here, /in is the reduced 
mass of the regular proton-electron system. In terms of 
the relative velocity between the two dark atoms, the 
cross section reads 
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FIG. 23. Momentum-transfer cross section for hydrogen 
156 (purple squares) and positronium (Ps) |157| (blue disks) 
self-scattering. We have rescaled the energy dependence of 
the H-H cross section according to Eq. (D12l to compare 
it to the positronium cross section. We also display the 
momentum-transfer cross section we use for dark-atom scat- 



tering (Eq. (D15l) 



where we used E = m^v 2 /A. In the limit where 



m p , Eq. (D16) agrees within a factor of order unity with 



the positronium (Ps) self-scattering cross section [157] 
(see Fig. 23 ). At very low velocities (v < 10~ 4 ), the above 
expression likely overestimates the momentum-transfer 
cross section. However, since very little energy is trans- 
fered in these low velocity collisions, we do not expect 
this to have a large impact on our results. 



2. H-p and H-e Cross Section 

When a dark proton collides with a neutral dark atom 
such that no atomic transitions are excited, two distinct 
processes must be taken into account. The dark proton 
can either collide elastically on the dark atom or it can 
capture the dark electron in a symmetric charge trans- 
fer Similarly, when a dark electron encounters a neutral 
atom at low energy, it can either scatter elastically or 
eject the atomic dark electron and be captured by the 
atomic nucleus. Since the final products of these two 
types of processes are indistinguishable, they both need 
to be taken into account to accurately calculate the dif- 
ferential elastic cross section |135j 
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where /d(#) is the amplitude for the direct elastic scatter- 
ing while / cx ($) is the amplitude for the charge-exchange 
process. Quantum interference between these two types 
of elastic scattering generally implies that a mt p,e 7^ 



H-p,e 



Much of the above discussion for the atom-atom case 
applies here and we can relate the dark proton-atom and 
electron-atom momentum-transfer cross section to that 
of standard proton-hydrogen and electron-hydrogen scat- 
tering as 
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where /i_D P , e = TO DW P)e /(m£> + TO p,e) and similarly for 
/^Hp.c- From the results presented in Ref. |135j . we 
can obtain an approximate expression for the proton- 
hydrogen momentum-transfer cross-section (obtained 
through a numerical fit) 
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again valid for 10 < Ecm/B^ < 10. In the case 
of electron-hydrogen scattering, the momentum-transfer 
cross section is smaller by almost two orders of magnitude 



<7^°{E C m « B H ) w 7na 2 . 
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The main reason for this large difference is that the 
electron-hydrogen differential cross section is strongly 
peaked in the forward direction, which implies small 
momentum transfer. One the other hand, the proton- 
hydrogen differential cross section has a significant back- 
ward peak which contributes substantially to the overall 
momentum-transfer cross section. Phcnomcnologically, 
we can take this into account by rescaling the cross sec- 
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Eq. D23 
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in contrast to the symmetric atom-atom case. 



FIG. 24. Momentum-transfer cross section for electron- 
hydrogen scattering obtained from direct computation |158l 
I159| (purple squares), from rescaling the proton-hydrogen 
cross section [156] (blue disks), and from rescaling the 
electron-positronium cross section [160] (yellow diamonds). 
We also display the momentum-transfer cross section we use 
for dark electron-atom scattering (Eq. ( |D23| |). 
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tions by the ratio of the masses of the two colliding bod- 
ies, that is 



or in terms of the relative velocity between the ions and 
the dark atom 



H-p 



2m r 



(D21) 



where we have used Eq. (D12) to appropriately rescale 



the collision energy. To obtain insight on the case m e ~ 
m p , we also consider the case of an electron elastically 
scattering off positronium |160j . In this case, we expect 
the low-energy cross section to be somewhat suppressed 
since the first Born approximation for the scattering am- 
plitude exactly vanishes. Existing computations support 
this fact and predict a roughly constant cross section at 
low energy. With this in mind, we conservatively neglect 
the small energy dependence of the H-p cross section at 
low energy and write 
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where we have used E = fiae^v /2. In Fig. 24 
we display the H-e momentum-transfer cross section, 
both from direct computation and from rescaling the p- 
H and electron-positronium cross sections according to 
Eq. (D21). We see that this rescaling is accurate up to 



a factor of order unity over the energy range of interest 
and thus provide the correct order of magnitude for the 
cross section. We also display in the figure our approxi- 
mate expression for the dark electron-atom cross section 
Eq. (lD23|). 
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